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This article presents the first measurement of the ratio of branching frac- 
tions B (A° N^ pTVf/j /B (A° A^TT^). Measurements in two control sam- 
ples using the same technique, B {B^ — > D+/i~i^j^) /B {B^ — > TJ+tt") and 
^ Z)*(2010)+/u"i/^) ^ Z)*(2010)+7r-), are also reported. The analysis uses 
data from an integrated luminosity of approximately 172 pb~^ of pp collisions at ^/s = 1.96 TeV, 
collected with the CDF II detector at the Fermilab Tevatron. The relative branching fractions are 
measured to be: 

B{K^,-^K,^i i/^) ^ ^ ^ +2.6 ^ q 3 

B {B° D+n-f^ 



B {B° D+TT-) 

B{B" ^ P- (2010) ji-i^^) 
B{BO -> D*(2010)+7r-) 



= 9.9 ± 1.0 (stat) ± 0.6 (syst) ± 0.4 (PDG) ± 0.5 (EBR) , 
= 16.5 ± 2.3 (stat) ± 0.6 (syst) ± 0.5 (PDG) ± 0.8 (EBR) . 



The uncertainties are from statistics (stat), internal systematics (syst), world averages of measure- 
ments published by the Particle Data Group or subsidiary measurements in this analysis (PDG), 
and unmeasured branching fractions estimated from theory (EBR), respectively. This article also 
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presents measurements of the branching fractions of four new A° semileptonic decays: A^J 
Ac (2595)+ n-u^, Ag ^ Ac (2625)+ ij-'u^, Ag ^ Ec{2i55)°TT+ fj.- u^, and AJJ ^ Ec(2455)++7r"/i'i^^, 
relative to the branching fraction of the A° A+/^~i?^ decay. Finally, the transverse-momentum 
distribution of AjJ baryons produced in pp collisions is measured and found to be significantly dif- 
ferent from that of mesons, which results in a modification in the production cross-section ratio 
a^o /(Jgo with respect to the CDF I measurement. 

PACS numbers: 12.39.Hg;13.20.-v;13.20.He;13.25.-k;13.25.Hw;13.30.- 

a;13.30.Ce;13.30.Eg;13.60.Rj;13.75.Cs;13.85.-t;14.20.-c;14.20.Mr;14.40.-n;14.40.Nd;23.20.Ra;25.75.Dw 
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Amplitudes for the weak decays of b hadrons are de- 
scribed by the product of Cabibbo-Kobayashi-Maskawa 
(CKM) matrix elements P, [3] and dynamical fac- 
tors. The CKM matrix elements represent the coupling 
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strength of the weak decays and are fundamental pa- 
rameters of the standard model of particle physics. In 
order to extract values of the CKM elements, knowledge 
of the dynamical factors is needed either from experi- 
ment or theory. Calculation of the dynamical factors, 
in the case of 6-hadron decays, is aided by heavy quark 
effective theory (HQET) HQET is an approx- 

imation relying on the large mass of the b quark (mf, w 
4000 MeV/c^), as compared with the Quantum Chro- 
modynamics (QCD) energy scale (Aqcd ~ 200 MeV), to 
imply a spin-independent interaction between the b quark 
and the light degrees of freedom. In baryon spectroscopy, 
the light degrees of freedom are in a relative spin-0 state 
for all A-type baryons; there is no spin-related interac- 
tion between the b quark and the light degrees of free- 
dom. Therefore, the subleading order corrections to the 
heavy quark limit are simpler than those mesons which 
contain a b quark {b mesons) '6|. Measurements of A°- 
baryon branching fractions may be compared with pre- 
dictions by HQET and test the calculation of dynami- 
cal factors to subleading order. However, in contrast to 
the b mesons, little is known about the A" baryon. At 
the time of writing this article, only five decay modes 
of the A° have been observed, with large uncertainties 
on their branching fraction measurements Q. On the 
theoretical side, combining measurements of the CKM 
matrix element, \Vcb\, and the world average of the A° 
lifetime 0, [1] , the branching fraction predicted by HQET 
for A° A+/z~P^ is 7.6% by Huang et al. M, and that 
for A° A+TT" is 0.54% by Leibovich et alP^. An in- 
dependent prediction of S(A" A+7r~) by Cheng using 
the non-relativistic quark model yields 0.50% [11.]. 

Presented here is the first measurement 
of the ratio of the A° branching fractions, 
B (AJJ A+/i~j>^) /B (Ag A+TT"). This measurement 
is based on data from an integrated luminosity of ap- 
proximately 172 pb~^ of pp collisions at y/s = 1.96 TeV, 
collected with the CDF II detector at Fermilab. Taking 
advantage of the relatively long lifetime of b hadrons 
(cT Ki 400 /xm), all 6-hadron decays described in this 
article are reconstructed from data satisfying an online 
event selection (trigger) which requires two charged 
tracks forming a vertex displaced from the location of pp 
collisions (two-track trigger). The A+ is reconstructed 
using the three-body decay A+ — > pK~Tr^, therefore 
both the — > A+7r~ and A+/x~i/^ decays 

result in four charged particles which are observable 
in the detector and have a similar topology (Fig. [T]). 
Since both decays have a similar topology and satisfy 
the same trigger, most systematic uncertainties from 
the detector, trigger, and reconstruction efficiencies 
cancel in the measurement of the ratio of branching 
fractions. Throughout this article, the inclusion of 
charge conjugate decays is implied. 

The ratio of the branching fraction of the A" exclusive 
semileptonic decay relative to that of the A" hadronic 
decay, Scxciscmi/^had, is extracted from the ratio of signal 
yields (A^cxciscmi/Miad) divided by the ratio of acceptance 



times efficiency (eexcisemi/chad): 

^ i^b ^ 1^ ^m) _ ^cxclsemi 
B{Al^A+Tl-) " Shad 

(-^exclscmi \ / / ^cxclsemi \ 
A'^had / / V fhad / 
I A^inclscmi A^scmibg | ^had ^-|^^ 
\ A^had / ^cxclsemi 

The analysis procedure can be summarized in four steps. 
First, the hadronic (A^7r~) and inclusive semileptonic 
(A+^^X) candidates are reconstructed. Second, the 
yields, A^had and Mnciscmi, are determined by fitting the 
mass distributions. Third, the contribution of back- 
grounds which produce a A^//~ in the final state is either 
measured or estimated and combined into A'somibg- The 
estimate of A'somibg requires a modification of the pro- 
duction cross-section ratio, crf^o/ago, with respect to the 
CDF I measurement [l^]. The dominant backgrounds 
which contribute to A'somibgi A" A^ (2595)^ 
Ag ^ A, (2625)+ /^-P^, AO ^ S,(2455)0^+/^-^^, and 
A° Sc(2455)++7r^^^i^^, have also been reconstructed 
in the data for the first time. Measurements of their 
branching fractions relative to the branching fraction of 
the A° — > A+^^t'^ decay will be used in the estimate 
of A'^gomibg- Fourth, the ratio of the products of detector 
acceptance and reconstruction efficiency, ehad/eexcisemi, is 
estimated from simulation. 

The analysis method described above is tested 
by performing the same measurements in de- 
cays which have a similar event topology. Specifi- 
cally, the following ratios of branching fractions are 
measured: B (S" ^ D+ii-D^,) jB {B° D+tt-) 

where — > K^tt'^tt^, and 

B D*{2010)+fi'Df,) IB D*(2010)+7r-) 

where D*(2010)+ -> D^t:+,D^ -f K-7r+. The results 
of the measurements are compared with previous 
results from the B factories to check the techniques 
used in this analysis. 

This article is structured as follows. Section |lT] de- 
scribes the relevant parts of the CDF II detector and trig- 
ger. Section [HT] details the event selections for the A+7r~ 
and A'^n~X samples. Section ITVl describes the simula- 
tions used in this analysis. Section |V] gives an account 
of the determination of the yields Miad and Mnciscmi- In 
Section IVIl A'^gomibg is estimated. Section IVIII includes 
measurements and estimates of the branching fractions 
of other A^ semileptonic decays which may contribute to 
A^scmibg, and an estimate of S(A° A+7r~) derived from 
a modification of the CDF I measurement of a/^o /a^a . 
Section IVIIII summarizes the systematic uncertainties. 
Section IIXI shows the measurements with the B'^ control 
samples using the same analysis technique. Section |X] 
compares the results of the A^ and -8° relative branch- 
ing fractions with the predictions from HQET and the 
world averages, respectively. Finally, Section lXll gives the 
conclusion. Unless stated otherwise, branching fractions. 
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fragmentation fractions, and lifetimes are obtained from 
the Particle Data Group world averages • The symbols 
"-ffc" and "TJfc" are used to generically denote hadrons 
containing charm and bottom quarks, "c hadrons" and "5 
hadrons" , respectively. The symbol "MC" , which stands 
for "Monte Carlo" , is used to generically denote simula- 
tion. 



II. THE CDF II DETECTOR AND TRIGGER 

The CDF II detector is a cylindrically symmetric ap- 
paratus described in detail elsewhere ^13;] . Only the parts 
of the detector relevant for this analysis are summarized 
here. The crucial features of the detector for this mea- 
surement are the tracking and muon systems. The track- 
ing system, which enables reconstruction of the trajecto- 
ries of charged particles, is contained in a superconduct- 
ing solenoid which gene rates a 1.4 Tesla magnetic field 
in the — z direction [14| . The 96-cm long silicon vertex 
detector (SVX II) [1^ consists of six equal subsections 
in z and five concentric layers of double-sided silicon sen- 
sors from r = 2.45 cm to r = 10.60 cm. The 310-cm long 
central outer tracker (COT) [l^, an open-cell wire drift 
chamber, consists of 96 sense wire layers from r = 40 cm 
to r = 137 cm which are grouped into alternating axial 
and ±2° stereo superlayers. The SVX II and COT pro- 
vide both r-0 and z position measurements in the pseu- 
dorapidity region oi\r]\ < 2 and |»7| < 1 respectively. 
The 452-cm long central muon detector (CMU) [l^j a set 
of drift chambers mounted outside of the central hadron 
calorimeter at r = 347 cm, contains 4 sense wire lay- 
ers which allow the formation of short track segments 
(stubs) and identify the muon candidates in the region 
of |77| < 0.6. 

The data for this analysis are collected with a three- 
level, two-track trigger which selects events with a dis- 
placed vertex. Consequently, data satisfying this trig- 
ger are rich in heavy flavor with a low background from 
the combination of random tracks (combinatorial back- 
ground). A schematic diagram of the event topol- 
ogy and trigger requirements is shown in Fig. [T] The 
strategy of the two-track trigger is as follows: at the 
first trigger level, the extremely fast tracker (XFT) [l^ 
finds two oppositely-charged tracks reconstructed in the 
COT, with a minimum transverse momentum (px) of 
2.04 GeV/c for each track. The scalar sum of the pr from 
the two tracks is required to exceed 5.5 GeV/c and the 
azimuthal angle between the two tracks (A0) to be less 
than 135°. At the second trigger level, the silicon vertex 
trigger (SVT) Q attaches hits measured with SVX II to 
the tracks found by XFT. The SVT reapplies the px re- 
quirements made at level 1 and further requires that each 
track has a transverse impact parameter (dp), measured 
at the point of closest approach (POCA) with respect to 
the beam line [2lj, in the range 120 /im-1000 /im. In 
addition, A(j) between the two trigger tracks is required 
to be in the range 2° - 90°. The intersection of the two 
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FIG. 1: An r-0 view of a Aj^ A+tt" (Ag ^ A-tfJ.~i^^) de- 
cay with a two-prong Aj decay vertex and a three-prong A^ 
decay vertex. In this case, the do of each pion (the pion and 
the muon) and the L^^^^^ of the two-pion vertex (pion-muon 
vertex) satisfy the trigger requirements. 



tracks forms a displaced vertex. Finally, the quantity, 
L"^^^^^, defined as the projection of the vector from the 
primary vertex to the displaced vertex onto the vector of 
the total momentum of the two tracks in the r-4> plane, 
must be larger than 200 /im. The level 1 and 2 triggers 
are implemented in hardware, while at the third level, a 
cluster of computers uses all detector information to per- 
form a full reconstruction of the event [2^ . In addition to 
reinforcing the same requirements as applied at level 2, 
level 3 requires the difference in z between the two tracks 
at the POCA to be less than 5 cm. The measurements 
presented in this article are based on an integrated lu- 
minosity of w 172 pb~^ collected between February 2002 
and September 2003, comprising w 152 million two-track 
trigger events. 



III. EVENT RECONSTRUCTION 

The final states A+7r~ and A^^~X, where A+ 
pK~TT^, are reconstructed in the data collected with the 
two-track trigger. The selection criteria for the hadronic 
and the inclusive semileptonic decay modes are kept as 
similar as possible, which reduces systematic uncertain- 
ties on the relative branching fractions. 

Both signal decays have a four-track topology. There- 
fore, events are selected which contain a minimum of four 
tracks, each with a minimum px of 0.5 GeV/c, do less 
than 5000 /im (measured with the SVX II), a minimum of 
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20 hits each in the COT axial and stereo layers f23j , and a 
minimum of three axial hits in the SVX II. Each track is 
also required to be in the fiducial region of the COT and 
to traverse all 96 wire layers. Making these requirements 
on each track ensures good quality of the track recon- 
struction and good momentum resolution. In addition, 
the maximum requirement on do suppresses background 
from daughters of Kg and A°, and from particles pro- 
duced by inelastic collisions of beam products with the 
detector material. 

The reconstruction begins by identifying the A+ can- 
didate. Only combinations of three tracks which satisfy 
the requirements described above are considered. Ev- 
ery combination must have two positively charged tracks 
and one negatively charged track. At least one of the 
three tracks must match a displaced track found by the 
SVT (SVT track [11]). The proton mass is assigned to 
the positively charged track of higher px, the pion mass 
to the track of lower pT, and the kaon mass to the nega- 
tively charged track. Assuming the proton track to be the 
higher pT track reduces the combinatorial background by 
« 50% while keeping « 90% of the A+ signal. A three- 
track kinematic fit determines the A+ decay vertex by 
varying the track parameters of the daughter particles 
simultaneously, within their uncertainties, so that the 
between the adjusted and the original track parameters is 
minimized. Only three-track candidates for which the fit 
converges and the invariant mass {MpK-w) is in the range 
2.18 ~ 2.38 GeV/c^ are considered further. 

Next, the selected A+ candidate is combined with an 
additional negatively charged track to form a A" can- 
didate. This fourth track must be matched to an SVT 
track. The combination is considered a A'^ semileptonic 
candidate and a muon mass is assumed for this track 
if the following two requirements are satisfied. First, a 
CMU muon stub must be present within 30 cm of the 
extrapolated track at the CMU radius (r = 347 cm) in 
the r-<j) view. Second, the matching between the track 
and the stub positions '25] is less than 9. Otherwise, the 
combination is a A^ hadronic candidate and a pion mass 
is assumed. Both the muon and the pion tracks from 
the A° decay must extrapolate to the fiducial region of 
the CMU. Making the same fiducial requirement for the 
hadronic and semileptonic modes ensures that the track- 
ing efficiency of both modes cancel in the ratio. 

Once all four A^-candidate tracks are found, the two 
tracks which have been matched to SVT tracks (one track 
from the A+ candidate, the other is the fourth track) 
must pass the two-track trigger requirements as described 
in Section [Til Then, a four-track kinematic fit is per- 
formed. This fit includes two constraints. First, the 
daughter tracks of the Ajj" must originate from a com- 
mon, tertiary vertex. Second, the trajectory of the A^ 
candidate must intersect with that of the remaining A^- 
candidate track, in three dimensions; this intersection 
is the decay vertex of the A|J candidate (defined as the 
secondary vertex). The secondary and tertiary vertices 
are determined in the four-track kinematic fit simulta- 



neously. These constraints improve the precision of the 
A+ decay vertex determination and the invariant mass 
of the A+ candidate is recalculated. After the kinematic 
fit, the values of MpK-w must be in the range: 2.269 - 
2.302 GeV/c^ {2a around the mean) for the hadronic can- 
didates and 2.18 - 2.38 GeV/c^ for the inclusive semilep- 
tonic candidates (see Fig. [2]). The wider A+ mass window 
for the semileptonic candidates allows for the Mpx-n spec- 
trum to be fit to extract the yield A'^incisemi- Also for the 
semileptonic decays, the four-track invariant mass Mj^^^ 
must be in the range of 3.7 - 5.64 GeV/c^, where the 
minimum requirement on Ma^^ reduces the background 
from other c-hadron and 5-hadron decays. See Section IVTl 
for more details. 

In order to reduce the combinatorial backgrounds fur- 
ther, the selection criteria on the following variables are 
optimized: px of the proton track, px of the fourth A^- 
candidate track [px (tt", /i~)], px of A+, px of the four- 
track system, of the A+ and the four-track kinematic 
fits, proper decay length (ct) of the A+ candidate, and 
(pseudo) proper decay length (ct*) of the A° candidate. 
The Xr,p is the r-(f) plane contribution to the returned 
by the kinematic fit. The ct is defined as: 



ct : 



(2) 



where is the projection of the vector from the sec- 
ondary to the tertiary vertex onto the momentum vector 
of A+ in the r-(j) plane, and is the world average of 
the A+ mass 171. The ct* has a similar definition: 



ct* 



Mao 



^'■^PT (4trks)' 



(3) 



where L^^ is the projection of the vector from the pri- 
mary to the secondary vertex onto the total momentum 
vector of the four tracks in the r-(p plane, px (4trks) is the 
transverse component of the total momentum of the four 
tracks, and M^o is the world average of the A° mass 0- 
Here, the primary vertex is estimated from the intersec- 
tion of the beam line and the trajectory of the A° candi- 
date. 

The optimization procedure maximizes the signal sig- 
nificance of the hadronic decays, S/y/S + B, where S is 
the number of A'^ — > A^tt" events in simulation multi- 
plied by a data-to-MC scaling factor and B is the number 
of background events estimated from the A+tt^ candi- 
dates in the data sideband. The data-to-MC scaling fac- 
tor for S is obtained by comparing the number of signal 
events in data and simulation with relaxed requirements. 
The background B is estimated by fitting the mass side- 
band region above the A^ signal peak with an exponential 
function and then extrapolating from the sideband region 
to the 3tT signal region around the peak. The optimized 
selection criteria are listed in Table U Figure [DJ a) shows 
the reconstructed A/acTt spectrum from the hadronic data 
and Fig. ^h) shows the reconstructed MpKn spectrum 
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TABLE I: Optimized requirements for reconstructing the 
A° A^TT^ and Aj/i^X decays. 



Pt (p) 


> 2 GeV/c 


PT (tt",/^") 


> 2 GeV/c 


PT (A+) 


> 5 GeV/c 


PT (4trks) 


> 6 GeV/c 


X?,(A+) 


< 14 


Xr0 (4trks) 


< 15 


rf(A+) 


> -70 /im 


cr (AO) 


> 250 



from the inclusive semileptonic data, both after applying 
the optimized selections. The most significant peaks in 
Fig. [5] represent the signals for each decay mode. In or- 
der to obtain the correct signal yields, a good modeling of 
the mass spectra, which includes a description of signal 
and background, is needed. The mass spectrum shapes 
of backgrounds from partially-reconstructed or misiden- 
tified 6-hadron decays are determined by fitting the mass 
distributions from simulation. The next section describes 
details of the simulations used in this analysis. 



IV. SIMULATION 

In order to determine the mass spectrum shapes close 
to the signal peaks in Fig. [2] and to estimate the ac- 
ceptance and efficiency of signal and background, both 
generator- level and full simulations are used. The 
generator-level simulation includes only the production 
and decay of b hadrons, and the analysis requirements are 
applied to quantities immediately after generation. The 
full simulation includes simulation of the CDF II detec- 
tor and trigger, and track reconstruction. It was found 
that the efficiency ratios, ehad/ecxciscmi, from a generator- 
level simulation and from a full simulation differ by only 
« 3%. The generator-level simulation is used to estimate 
the quantities which are found to be small or already 
have large uncertainties from other sources [2^ : the size 
of the background contribution where the A^ and the 
fi~ originate from two different heavy-flavor hadrons pro- 
duced by the fragmentation of bb or cc pairs (termed bb/cc 
background) , part of the AJJ systematic uncertainties (the 
semileptonic decay model and lifetime of A°, A° and A+ 
polarizations, and Ajj" Dalitz structure), and modification 
of the CDF I (Jf^a/dgo result. Therefore, this 3% differ- 
ence has a negligible effect on the final measurement. The 



following subsections describe the key components of the 
simulations used in this analysis. 

A. Production and decay of b hadrons 

Two different programs are used to simulate 6-hadron 
production: pythia version 6.2 [23|, which simulates all 
of the strong interaction processes which are involved in 
6-hadron production, and bgenerator ^28^, which gen- 
erates a single h hadron in the event. Since pythia simu- 
lates all of the products of the pp collision, it is computa- 
tionally intensive to produce a given final state. There- 
fore, PYTHIA is used to estimate only the bb/cc back- 
grounds in the inclusive semileptonic data f Appendix [Cjl . 
The PYTHIA generator simulates physics processes using 
leading-order matrix elements, supplemented by initial 
and final state radiation. The program also includes 
hadronization of the quarks and gluons in the final state 
and the beam remnants left when a parton undergoes 
high-momentum scattering. The bgenerator program 
is very efficient at producing a large sample of a specific 
6-hadron under well-defined kinematic conditions. It is 
used to determine the acceptance and efficiency for signal 
and other backgrounds and to model the mass spectra. 
In the BGENERATOR program, a single h hadron is gen- 
erated using the measured px spectra of 6 hadrons as 
inputs. The and pt spectra are derived from the 
fully-reconstructed A^ A+7r~ and i?" — > ZJ+tt" de- 
cays in the two-track trigger data, after correcting for 
acceptance and efficiency. 

After the event generation, the decays of the b and c 
hadrons and their daughters are simulated using the evt- 
GEN package [2§|. For all other particles in the event, this 
is done by the pythia program. The evtgen program 
uses the dynamics from a full matrix-element calculation 
and is tuned to measurements, mainly results from exper- 
iments at the T (45') resonance [13, Isil, IH, [11] where the 
decay models for the and the B~ have been demon- 
strated to match data. As a full theoretical model for A° 
semileptonic decays is not yet implemented in evtgen, 
a flat phase space (termed PHSP) simulation is used for 
A° decays. A correction is applied after generation to 
account for the proper A^ semileptonic decay dynamics. 
Details of this correction arc given in Section IIV CI 



B. Detector simulation and comparison of 
kinematic distributions 

After an event has been simulated at the generator 
level, it is processed with a full simulation of the CDF 
II detector and trigger. The geometry and response of 
the active and passive detector components are simulated 
using the GEANT software package [13]. The events are 
then processed with a two-track trigger decision program 
and reconstructed using the same executable as that used 
to reconstruct the data. The resulting events have the 
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FIG. 2: The reconstructed invariant mass spectra in data after applying all selection criteria: (a) the Ma^h spectrum of the 
hadronic candidates; (b) the MpKir spectrum of the inclusive semileptonic Kt candidates. 



same structure and format as the data and are analyzed 
in the framework described in Section IIIII 

Distributions of kinematic variables from the full sim- 
ulation with BGENERATOR input are compared with 
the same distributions from data. In order to com- 
pare the data and the simulation, the data distributions 
are background-subtracted. The agreement between the 
data and the simulation is quantified by a comparison 
probability and the ratio of spectra produced from the 
data and the simulation. All relevant distributions agree 
satisfactorily. Figure [3] shows good agreement between 
the data and the simulation in the pr (A") and pT (5") 
spectra. The pT of the h hadron is the most important 
kinematic variable in this analysis because the 6-hadron 
momentum is distributed among three particles in the 
final state for the exclusive semileptonic decay and be- 
tween two particles for the hadronic decay. 

C. Acceptance and efficiency scale factors 

In order to obtain accurate estimates of the acceptance 
and efhciency, several scale factors are applied to the 
number of events selected in simulation and their val- 
ues are listed in Table HTl As mentioned earlier, evtgen 
contains only a phase space simulation of semileptonic A" 
decays. In order to estimate the effect of decay models 
on the signal acceptance, a weighting of the flat phase 
space distribution according to a form factor model from 
Huang et al. @ for the hadronic current of the A|J to 
A+ transition, and a V — A model for the leptonic cur- 
rent, is performed at the generator level. The ratio of 
the generator-level acceptance after weighting relative to 



that before weighting, Cmodci, is found to be 0.994±0.025. 
Since this ratio is consistent with unity, the PHSP full 
simulation samples are used throughout the Aj] analy- 
sis. The correction factor, Cmodci, which accounts for 
the Al semileptonic decay dynamics, is applied to the ef- 
ficiencies for semileptonic decays. The shape of the Ma^i_i 
distribution is sensitive to the decay dynamics and may 
be used to cross-check the form factor and V — A mod- 
els (termed FF). Figure |4] shows the reconstructed Ma^^ 
distributions from the data and from the PHSP full sim- 
ulation, before and after multiplying the MC histogram 
with the bin-by-bin ratios which are derived from the 
same generator-level simulation samples for Cmodci [ssj . 
The corrected distribution has a significantly improved 
agreement with the data, which confirms the procedure 
for deriving Cmodci- 

In addition, the CMU muon reconstruction efficiency 
is found to be over-estimated in the full simulation; the 
resulting scale factor, Ccmtt, is measured using a sam- 
ple of J/ip decays The dependence of the 
XFT trigger efficiency on the particle species and pt is 
not included in the full simulation. Using a pure pro- 
ton sample from the A° —> p~^tt~ decays, and pure kaon 
and pion samples from the J* (2010)+ — > D'^tt^ decays 
where D° K-'k+ [H, [13, [ll], the data-to-MC scal- 
ing factors, Cp, Ck, and C^, are derived and applied to 
the track which passes the trigger requirements in the 
reconstruction program. 

With a reliable simulation for the modeling of mass 
spectrum shapes, the numbers of signal events can be 
determined by fitting the invariant mass spectra in Fig. [2] 
as described in the following section. 
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FIG. 3: Comparison of reconstructed pr spectra of 6-hadrons between the data and the fuU simulation; (a) A° and (b) S". 
The top figures show the pr distributions while the bottom figures show the ratio of data to MC. The curves in the bottom 
figures are the result of an exponential fit of the ratio. 



V. DETERMINATION OF THE SIGNAL 
YIELDS 

The numbers of hadronic events (-/Vhad) and inclusive 
semileptonic events (Mnciscmi) in Eq. ([1]) are extracted 
by fitting the Ma^^ and MpK-n spectra in data, respec- 
tively. The fit to the mass spectra is performed using 
an unbinned, extended likelihood technique [39], where 
the fit parameters are adjusted to minimize the negative 
log-likelihood (— ln>C). The general unbinned, extended 
log-likelihood is expressed as: 

\nC = ^ In [A^sig^ (m,) -f A^bgS (to,)] 

i 

-iVsig-A^bg + ^lnC,, (4) 

3 

where i represents the i^^ candidate and to represents the 
reconstructed mass AfA^Tr or MpK-n- The numbers of sig- 



nal and background events are denoted as A'sig and iVbg; 
S {rrii) and B {rrii) are the normalized functions which 
describe the shapes of signal and background mass spec- 
tra, respectively. Each Cj is a Gaussian constraint on a 
specific fit parameter Xj-. 

Cj = g {xj,fij,aj) = e ^\ J , (5) 

V ZTTCTj 

where the parameter Xj has a central value of fij and 
an uncertainty of aj. Because the data sample size is 
not large enough to determine these parameters accu- 
rately from the fit, they are constrained to values which 
are estimated from independent measurements and the 
full simulation. Definitions of the constrained parame- 
ters (xj) are given in Section IVAl 

Correct modeling of the mass spectra is crucial in 
the determination of iVhad and -/Vincisemi- Two types 
of background appear in each mass window of inter- 
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TABLE II: Acceptance and efficiency scale factors applied to 
the number of events selected in simulation. The pr is the 
transverse momentum (in GeV/c) of the track which passes 
the trigger requirements. The uncertainty on Cp is obtained 
by taking the difference between the pr-dependent formula 
below and a constant from an average over the A° sample, 
0.905. The uncertainties on Ck and Ctt are below 0.5% and 
have negligible effect on the final relative branching fractions. 



Scale Factor 



Value 



C*modcl 
CCMU 



p 



Ck 



0.994± 0.025 
0.986 ±0.003 
1.06 

0.969- 
1.002 - 



1J3 I 3.2 _ 2.2 

PT pfT 

0.094 

PT 

0.067 



PT 



est. The first is combinatorial background. Combina- 
tions of four random tracks contribute to tliis background 
in both the hadronic and semileptonic modes. Com- 
binations of a real c hadron and a random track con- 
tribute only in the hadronic mode. The mass spectrum of 
the combinatorial background is determined using data 
sidebands. The second background is misidentified or 
partially-reconstructed decays of b hadrons. Their mass 
spectrum shapes are determined using the simulations 
as described in Section IIVI The dominant contributing 



decays are identified with a generator-level simulation of 
inclusive 5-hadron decays, and are categorized accord- 
ing to their mass spectrum shapes. Decay modes with 
similar shapes are generated together using a full sim- 
ulation, with the number of generated events for each 
decay mode proportional to the fragmentation fraction 
times the branching ratio, and are parameterized by a 
single function. The functional form for each combined 
background spectrum is determined empirically to match 
the shape of simulated mass distribution. The parameter 
values of each function are obtained by fitting the sim- 
ulated spectrum. When fitting data, the values of the 
shape parameters are fixed while the normalization is a 
free parameter. 



A. The Ma^tt spectrum for the AjJ A+tt yield 

Figure [3a) shows the Ma^tt spectrum with the fit re- 
sult superimposed. The likelihood fit is performed in 
the mass window Ma^tt = 4.6 — 7.0 GeV/c^, whereas 
Fig. ini^a) shows a more restricted mass range near the 
signal peak. The — A+tt" yield returned by the fit 
is 179 ± 19. The signal peak at A/a^^ w 5.6 GeV/c^ 
is described by a Gaussian function. The width of the 
Gaussian is constrained in the fit to reduce the uncer- 
tainty on iVsig. The constrained width is the product of 
a data/MC scale factor and the width of the Ma^-w distri- 



X cr^ ^, which 



bution in the full simulation, {cr'^^'' / a^^) 
is 0.0231 ±0.0012 GeV/c2. The scale factor, a^^^^^/a^^, 
is obtained by comparing the width of the invariant mass 
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distribution in data with that in the simulated events, 
using the D+7r~ decay which has a similar topol- 

ogy and a larger data sample size. The combinatorial 
background is parameterized by an exponential (light- 
gray filled region), where the exponential slope is deter- 
mined by the A^tt" candidates in the mass region above 
5.7 GeV/c^. The functions which describe the mass spec- 
tra of backgrounds from the misidentified or partially- 
reconstructed 6-hadron decays are determined from the 
simulated mass distributions. 

Details of the background from the misidentified or 
partially-reconstructed 6-hadron decays follow. The dou- 
bly Cabibbo-suppressed decays A° A^K^, with a 
pion mass mistakenly assumed for the kaon, are indi- 
cated by the black filled region. The ratio of the num- 
ber of doubly Cabibbo-suppressed decays relative to that 
of the signal mode, N\^k /N^Aa-n, is fixed to 8% in the 
fit. This value is obtained from the world average of 
measurements in the B'^ modes. Fully reconstructed b- 
meson decays with misidentified daughters produce a dis- 
tinct peak at Ma^^ « 5.5 GeV/c^ (wavy-line region). 
The B° D^TT~ decays, where — > K^tt'^tt^ and 
one of the pions is reconstructed as a proton, contribute 
about 50% to this background. The background from the 
remaining partially-reconstructed 5-meson decays has a 
monotonically falling mass distribution (dark-gray filled 
region) dominated by i?'^ — > D^^^v^, — > D^p^ 
where p~ tt^tt^, and B^ — > D*(2010)+7r^ where 
Z3*(2010)+ Z3+7r° and the 7r"'s are not reconstructed 
in the event. The remaining decays also have a 
falling mass spectrum (hatched region) dominated by 
Al — > A'^iDi and A|J — > A+p" where p~ n°n~. 

B. The MpKiT spectrum for the A'^p~X yield 

Figure (HKb) shows the MpKTv spectrum for events with 
muons, with the fit result superimposed. The inclusive 
A+p-X yield returned by the fit is 1237 ± 97. The fit 
for the MpKiT spectrum is less complex than that for the 
-^AcTT spectrum described above. Note that the signal 
peak includes the backgrounds which also contain A+p^ 
in the final state (see Section IVI[) . The signal peak at 
MpK-K ~ 2.3 GeV/c^ is modeled by a Gaussian func- 
tion. Background from the 6-hadron semileptonic decays 
with a c-hadron daughter misidentified as a A+, such 
as i?'^ D'^p~Df^ where _fC~7r+7r+ and one of 

the pions is assigned the proton mass, does not produce 
a peak or distinctive structure and is inseparable from 
the combinatorial background. These two backgrounds 
are combined and modeled by a second-order polynomial 
(light-gray filled region) . 

C. Summary 

Table IIIII summarizes the A^ hadronic and inclusive 
semileptonic yields and the probability of correspond- 



TABLE III: Observed number of events in each AjJ decay 
mode determined from the unbinned, extended likelihood fit, 
X^/NDF, and the corresponding probability computed to in- 
dicate quality of the fit. 



Mode Yield xVNDF Prob (%) 
Aj]^A+7r- 179 ± 19 123/111 20.7 
A+^-A: 1237 ± 97 48/38 13.0 



ing fits. Each model describes the data well, as indicated 
by the probability. In order to obtain the number 
of exclusive semileptonic signal events A'excisemi, the con- 
tributions from backgrounds which also produce a A+ 
and a p~ in the final state, A^somibg: must be subtracted 
from A^inciscmi- Section lVTl describes the estimation of the 
composition of the A|J inclusive semileptonic data sample. 
Section IVIII details observations of four new A^ semilep- 
tonic decays and the estimates of A^ semileptonic and 
hadronic branching ratios which are required to deter- 
mine the sample composition in Section [VII 



VI. COMPOSITION OF THE INCLUSIVE 
SEMILEPTONIC DATA 

The B factories [s^, [3l|, [13, produce b hadrons in 
e~^e~ interactions where the beam energy may be used 
as a constraint when reconstructing events. This feature 
is particularly helpful for reconstructing semileptonic de- 
cays where a neutrino is missing. At the Tevatron, b 
hadrons are produced by the interactions between quarks 
and gluons with a broad parton momentum spectrum. 
Therefore, beam energy constraints are not available to 
aid 6-hadron reconstruction. Backgrounds which contain 
a A^, a p~ , and other particles in the final state can not 
be separated easily from the exclusive semileptonic sig- 
nal, A|^ A'^ p~Dp^, and will contribute to the inclusive 
A^ p~ events observed in data. These backgrounds arise 
from three sources: 

1. false muon: a A+ and a hadron track (p, or 
7r~) misidentified as a p~ . 

2. bb/cc: a A+ from the decay of a heavy-flavor 
hadron Hi, {He) and a /i^ from the decay of the 
other heavy-flavor hadron H\j {He), where the two 
hadrons are produced by the fragmentation of 66 
(cc) pairs. 

3. feed-in: decays of a single 6 hadron into a A+, a 
p~ , and particles not reconstructed in data. 

The goal is to measure the branching fraction of the ex- 
clusive semileptonic decay relative to that of the hadronic 
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FIG. 5: Results (curve) of the unbinned, extended likelihood fits for determining the numbers of AjJ candidates: (a) hadronic 
and (b) inclusive semileptonic. The filled histograms indicate various backgrounds. 



decay. The backgrounds listed above must be sub- 
tracted from the observed number of inclusive semilep- 
tonic events in data. Equation ([T]) is then re-written as 
follows: 

^cxclscmi 

Shad 

/ A^inclsomi - A^falsc^ - Nf^b,cc ~ ^^food \ Chad , 

M ■^^> 

\ J*had / ^cxclscmi 

The number of false- muon events (A'faiso/^) is obtained 
from data containing a A+ and a hadron track satisfy- 
ing reconstruction requirements, with the hadron track 
weighted by an appropriate muon-misidentification prob- 
ability. The contributions from the bb/cc (A^^f, ce) ^"^^ 
the feed-in backgrounds (Afoed) are estimated using both 
data and simulation. Instead of the absolute amount, the 
ratios Nf^i^^^/ N\^i^ and A'^fccd/Miad are estimated. Esti- 
mating the ratios instead of the absolute amount has one 
advantage: the majority of the background events are 
decays of A|J, so knowledge of the A^ production cross- 
section is not necessary. The quantities A^tb cc/-Ahad and 
A^feed/A'had are determined from the ratios of the prod- 
ucts of branching fractions and efficiencies (times produc- 
tion cross-section for non-Aj] background). The normal- 
ization procedure requires measurements or estimates of 
the branching fractions for the A" A+7r~ decay and 
for several semileptonic decays which may contribute to 
the backgrounds; details of these measurements and es- 
timates are found in Section IVlII The ratio Nf^i Ni^sd 
has been estimated to be very small and contributes 
< 1% to the A° A+^~j?p signal. More information 
on bb and cc backgrounds may be found in Appendix [C] 



The following sections describe the estimation of A'faiso/i 
and A^feod- 



A. False muons 

One type of semileptonic background is due to the 
pairing of a Ajf with a proton, a kaon, or a pion which 
is misidentified as a muon. A hadron is misidentified 
as a muon when it passes through the calorimeter into 
the muon detector, or when it decays into a muon in 
fiight. The probabilities for a proton, kaon, or pion to 
be misidentified as a muon {Vp, Vk, and T'^, respec- 
tively) are measured using a pure proton sample from 
the A*' — !• p+TT^ decays, and pure K and tt samples from 
the D*{2QIQ)+ -> D°t:+ decays where Z?" ^ K-'k+ 
The muon-misidentification probability is defined as the 
fraction of the CMU-fiducial and SVT-matched hadron 
tracks which satisfy the muon identification requirement 
(a track associated with hits in the CMU and with a 
matching less than 9). Figure [5] shows the Vp (mea- 
sured in twelve pt bins) and Vt^ Vk (measured in sixteen 
Pt bins) for positively and negatively charged tracks, sep- 
arately. A difference is observed between Vk+ and Vk- 
in the low pT region, which is not seen for protons and 
pions. The larger hadronic cross-section for the K^p 
scattering relative to that for the K^p scattering results 
in a lower rate of K~ being misidentified as muons pass- 
ing through the calorimeter. 

The contribution of the false-muon background to the 
^ A^A^^^M signal, A^faiscMJ is obtained by weighting 
data containing a A+ and a hadron track (h~), with the 
muon-misidentification probability (Pavg) as a function 
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of the momentum of h~ . This hadron track must ex- 
trapolate to the fiducial region of the CMU and fail the 
muon identification requirements in order to remove real 
muons. The other selection criteria for the A+ft,~ sample 



arc the same as those for the Aj] 



A+/i i^f^ reconstruc- 



tion. The TVfaisc^ is then extracted from a fit of the 
MpKTT distribution produced from the weighted A+/i^ 
sample. Figure [7] shows the result of the fit- 

Since no particle identification requirement is applied, 
whether h~ is a proton, a kaon, or a pion can not be de- 
termined from data. The muon-misidentification proba- 
bility, Pavg, is, therefore, an average ofVp, Vk, and Vtt 
weighted by Fp, Fx, and F^^ (the fractions of p, K, tt in 



FpVp + FkVk+F^V^. 



(7) 



In order to determine Fp, F^, and F^, physics processes 
which produce these hadrons must be understood. The 
principal sources of these hadrons after analysis require- 
ments are the decays Hi, ~* A^h~X, where h~ is a p, 
, or TT^ misidentified as a muon and X could be noth- 
A+7r~) or any other particles which are 



mg {e.g. 



not reconstructed {e.g.^ B~ — > K^p^~Vp). Other sources 
include: fragmentation of a primarily produced quark 
or gluon, inelastic collisions of secondary particles with 
the detector material, and decays of c hadrons. Hadrons 
which are not from &-hadron decays are suppressed by re- 
quiring that the transverse impact parameter (do) of the 
muon candidate is in the range 120 /Ltm-1000 /xm, and 
that the A^ and the muon candidates form a vertex sig- 
nificantly displaced from the beam line (see Section ITlT|) . 
In addition, the pythia simulation indicates that the 
background where a false muon and a A^ signal origi- 
nate from decays of two different 6 or c hadrons is less 
than 0.1% of the inclusive semileptonic signal and can be 
ignored. Therefore 
the Hb — > A+/i~X full simulation 

Table HVl shows values obtained for Fp, Fk, F^-, and 
-/Vfaiscp- The uncertainty on A'faiso/j includes: the sta- 
tistical uncertainty from the fit, the uncertainties on 
Fp, Fr-, and F^r, and the uncertainties on the measured 
muon-misidentification probabilities. The A^faisc/x is ap- 
proximately 3.2% of the number of A+/x^X events. 



B. Feed-in backgrounds 



Fp, Fk, and F^r are obtained from 



O b - 



6 
4 
2 


-2F- 



8 : 
6 : 
4 : 

2 \ 

-2 - 



P 
P 



■I, 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




||' m I'|'I' m I' m |' m I' m I'|| 



■71+ 
■71" 







( 


i 







5 10 15 

PT(GeV/c) 



FIG. 6: The probability for a proton, kaon or pion to be 
misidentified as a muon as a function of pr ^40.]. Note that 
for the measurements with negative values, a zero muon- 
misidentification probability is used to weight data. 



The feed-in backgrounds to the A° — > A+^ signal 
fall into three categories: 



1 Armcson 
^- ^:fcod 



Baryonic, semileptonic decays of 
W IB~ l&l^ which decay into A+, an anti-nucleon 
and leptons (e.i?., B~ K'^p^~Vp). 



T-b— baryon 

^0 ^A+KV-i^,.). 



2. iVf"„i''"''"": Semileptonic decays of other b baryons 



3. iV 



othcrA!^ ^ 
feed 



Other semileptonic decays of AP, which 



include either additional particles {e.g., A^ 
A^7r''"7r~/i~j>p) or a higher mass c baryon with 
subsequent decay into the A+ signal {e.g., A" 
A, (2595)+ ^i-D^, A, (2595)+ -> A+7), 



and the ratio iVfecd/A'^had is expressed as 

A^feed 



^^feed 



T\T-b— baryon . ^rOtherA^ 
^^tood + ^^feed 



^had 



(8) 
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FIG. 7: Distribution of MpK-w produced from data with a 
hadron track {h~) and a A J candidate in the final state, 
after weighting the hadron track with an average muon- 
misidentification probability (Pavg)- The curve indicates the 
result of the St- 



TABLE IV: The fractions of p, K, and tt in the h~ {Fp, Fk, 
and Ftt), the estimated number of false- muon events to the 
Aft fi~i/fj, signal, and for comparison, the number of the 

inclusive semileptonic events in data. 





0.24 ±0.16 


Fk 


0.05 ±0.08 


F^ 


0.71 ±0.16 




40 ±9 


A^iiiclscmi 


1237 ±97 



The TVi""!"" and N°~^'''^°'' have been estimated to be 



•feed '^^^^ ^*fcod 

very small and contribute < 1% to the A" 



A^ 



/i Vf, sig- 



nal. Details of these estimates are found in Appendices El 
and[B] The contributions from other Aj] semileptonic de- 
cays are estimated below. 

rOthcrA!! 



The ratio N, 



otherA 
feed 



feed "/A'^had is given by: 



ad 



B (AO ^ A^ 



(9) 



where Bi and are the branching fraction of A|] semilep- 
tonic decay mode i and the efficiency of partially re- 
constructing the decay i as the semileptonic signal, re- 
spectively. The estimate of Nf^^^'^^^ starts by identify- 



ing the dominant background decay modes that enter 
Eq. dH). The observation of spin-1/2 A^ (2595)+ and 
spin-3/2 Ac (2625)+ [ill, IH indicates the existence of 



cays. 



Ac (2595)+ /I- 



and Ag 



Ac(2625)+/z-p^ de- 



In addition, the following decays may contribute 



Ag 
Ag 
Ag 
K 



S]e(2455)0^+M-i>p, 
S]e(2455)+^0M-i>^, 
S]e(2455)++7r-/i~j>^, 
A+/o(980)Ai-p^, 
A+TT+TT^/i^i?^ (non- resonant), 
A+TrOTT^/i^P^ (non-resonant). 



The decay in the tau channel, A° A+r^j/,- where 
^^v^Vr, also makes a small contribution. Equation ^ 
requires knowledge of the branching fractions of A^ 
A+7r~ and these background decays. In order to reduce 
systematic uncertainties from theoretical predictions, the 
dominant background decays, A° Ac (2595)+ ^"i?^, 
Ac(2625)+/i-p^, Ag -> Sc(2455)"7r+/i-j>„, and 



Ag 
Ag 



have been reconstructed in 



(2455)++7r-^''-i^^! 
the data. Measurements of their branching fractions rela- 
tive to the branching fraction of the AJJ A+/i~P^ decay 
and estimates of the branching fractions of AjJ — > A+tt^ 
and the other A" semileptonic decays are found in Sec- 
tion IVIII Once the list of background decay modes is 
established and their branching fractions are estimated, 
the acceptances and efficiencies of these backgrounds rel- 
ative to that of the hadronic mode (ci/e^o^A+Tr- ) 
obtained from the full simulation as described in Sec- 
tion |TVl Figure [8] shows that a minimum requirement on 
-^^AcM of 3.7 GeV/c^ reduces the backgrounds from other 
A^ semileptonic decays which have more particles in the 
final state. 

Table fVl summarizes the feed-in backgrounds from the 
A" semileptonic decays discussed above and lists the 
hadronic and inclusive semileptonic yields observed in 
data. The two leading backgrounds after all selections 
are Ag ^ A^ (2595)+ m"!?^ and A^ ^ Ae (2625)+ i?^. 
The total contribution from feed-in backgrounds has been 
estimated to be 24.0% of the number of A+/i^X events. 



C. Summary 

Table IVII lists the values of all the background vari- 
ables which enter Eq. ([6|) and summarizes the composi- 
tion of the inclusive A^/i^ sample. The dominant signal 
contamination is from the feed-in background. The sec- 
ond largest background arises from false muons. The 
smallest background source is bb/cc. The estimate of 

-^focd°'^^''/-^had requires knowledge of the branching frac- 
tion of each feed-in decay and also the hadronic decay 
A|J — ^ A+7r~. The next section details the measurements 
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TABLE V: Feed-in backgrounds to A^ A+ fi from other A° semileptonic decays. The "*" indicates decays which have 
been reconstructed for this measurement and seen in data for the first time (Section IVII A[) . The second column lists the 
estimated branching fractions from Section IVIII Numbers in parentheses are estimated uncertainties [isl . [i^ . The third 
column lists ei/e^o^A+ir- ^i^^i statistical uncertainty. All efficiencies are determined from the full simulation as described 
in Section HVl The fourth and the fifth columns list the normalization for each background relative to the hadronic and the 
exclusive semileptonic signals, respectively. The last column lists the number of events for each background after multiplying 

7Vj°*j^°'^*''/7Viiad) by Aiiad, and the uncertainty includes only the statistical uncertainty on A^had- Note that while the numbers 



listed in the fourth column are used in the final measurement, the last two columns are shown only for a comparison with the 
A|J Afjj.~i>^ signal. 



Mode B{%) (^ife^j iV..e„. 

A" ^ A+TT- 0.36+[p-^gj 1.000 - - 179 ±19 

A+/i-A: - - - 1237 ±97 

Ag^A+^~t^^ 7.3±(1.4) 0.303±0.004 6.118 1.000 

* Ag ^ Ac (2595)+^-;^^ 0.9±(0.4) 0.198 ±0.003 0.503 0.082 90 ±10 

Sc(2455)++7r- 24 ±7 

^ A+Tr+ 100 ±(5) 

Sc(2455)°7r+ 24 ±7 

^ A+n" 100 ±(5) 

Sc(2455)+7rO 24±(1.2) 

^ A+tt" 100 ±(5) 

-^A+TT+vr- 18 ±10 

-^A+vrOvrO 9 ±(0.45) 

-^A+7 1±(0.05) 



* Ag^ A,(2625)+/i-i^^ 


1.5 ±(0.6) 0.192 ±0.003 0.815 


0.133 


146±15 


^ A+vr+vr^ 


66 ±(3.3) 






^ A+vrOvrO 


33±(1.7) 








1±(0.05) 











0.39 ±(0.23) 0.082 ±0.004 


0.089 


0.015 


16± 


2 




^A+n- 


100 ±(5) 










Ar 


Sc(2455)+7r°/i-p^ 


0.39 ±(0.23) 0.073 ±0.004 


0.080 


0.013 


14± 


2 




A+TfO 


100 ±(5) 










* AO- 


^ E,(2455)++7r-/i-p^ 


0.39 ±(0.23) 0.077 ±0.004 


0.084 


0.014 


15± 


2 




^A+n+ 


100 ±(5) 












-^A+T-Vr 


2.0 ±(2.0) 0.041 ±0.003 


0.040 


0.006 


7± 


1 






17.36±0.05 










Ar 


A+f 0(980) fi-D^ 


0.00 ±(0.32) 0.023 ±0.002 


0.000 


0.000 







Ar 


-> A+vr+TT^/i^Pp 


0.00 ±(0.64) 0.032 ±0.002 


0.000 


0.000 







Ar 


^ A+7r"7rV-i^M 


0.00 ±(0.32) 0.033 ±0.002 


0.000 


0.000 








and assumptions used to estimate these branching frac- VII. OBSERVATIONS OF FOUR NEW A^ 

tions. SEMILEPTONIC DECAYS AND ESTIMATES OF 

Aj; SEMILEPTONIC AND HADRONIC 
BRANCHING FRACTIONS 



The size of the background contribution from the feed- 
in of other semileptonic decays of A°, A^f°o'd''^'', is nor- 
malized to the observed hadronic signal yield in data, 
with corrections for the relative acceptance times effi- 
ciency for each decay mode [see Eq. ©J. This proce- 
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TABLE VI: The values of background variables in Eq. ([6]) and the composition of the A^/i X sample. Uncertainties on the 
bb/cc and feed- in backgrounds to the — > decay are statistical only. The values of A^'had, A'^indscmi, and ehad/eoxcisomi 

in Eq. (O are 179 ± 19, 1237 ± 97, and 0.303 ± 0.004, respectively. 



Source TV iV/iVhad ^Mnclscmi (%) 

Signal - - 72.5 

False muon 40 ± 9 - 3.2 

hh/cc ~ 0.017 ±0.009 0.3 

Fccd-in - 1.660 ± 0.018 24.0 
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FIG. 8: The fully-simulated Ma^/j distributions of the signal 
and backgrounds from other semileptonic decays of A°, on a 
log scale. The arrow indicates the minimum analysis require- 
ment at 3.7 GeV /(?. Note that the backgrounds tend to be 
at the lower mass region. 



dure requires estimates of the branching fractions for 
each background decay and also for the hadronic sig- 
nal. In order to reduce systematic uncertainties from 
these branching fractions, resonant semileptonic de- 
cays expected to contribute to the A+/x~ sample have 
been reconstructed in data. These reconstructed decays 
are then used to estimate the branching fractions of non- 
resonant semileptonic decays vi^ith a constraint of the 
viforld average of S (A° — ^ K'^l^Di anything) and an esti- 
mate of ;S(A° A+/x~P^). Note that the estimate of the 



A+7 



and AO 



A+/i branching fractions ap- 



pear in the estimate of the backgrounds (Section [VTB]). 
It must be pointed out here that from Table IVTl the to- 
tal contribution of feed-in background is 24.0% and con- 
tributes at most this amount to the total uncertainty for 
this measurement. Furthermore, the ratio of the esti- 



mated B (Ag A+pi-Dfj) to B (Ag A+TT-) need not 
be the same as the final measured result. 

Section IVII Al first presents measurements of the 

branching fractions of four new A^ semileptonic de- 



cays 



relative to that of the 



Ec(2455)°7r+^-p„,, and A? 



decay: 



A^ 

Ac (2625)'^ 
f,, c..^ ..J, ^ Sc(2455)++7r-/i-z>^, and 
then describes the estimation of B (A° A+zi^j?^) and 
branching fractions of several non-resonant A^ semilep- 
tonic decays. Section lVII BI shows that how the px distri- 
bution of the A^ baryon produced from pp collisions is sig- 
nificantly different from that of the meson and gives 
the corresponding modification to the ratio, cr^o/c^o, 
with respect to the CDF I measurement [l^l- The ra- 
tio ctao/ctb" is then used to estimate B (A" — > A+tt^). 



A. Observations of four new A° semileptonic 
decays and estimates of the Aj,' semileptonic 
branching fractions 

The following AJJ semileptonic decays are considered in 
the estimate of Affo!!d Section IVIBI 



Ag 
K 

K 
K 
K 
K 
K 
K 

Ag 



Ac (2595)+ 

Ac {2625)+ fi- 9^, 

S]c(2455)0^+M~i>;„ 

Sc(2455)+^0/^-z>^, 

Sc(2455)++7r-M"j>^, 

A+/o(980)Ai-p^, 

A+TT+TT^/i^P^ (non- resonant), 

A+Tr^TrO/i^p^ (non-resonant), 

A + T-D-r. 



Among the nine decays above, none have been observed 
previously [i^ ; only the branching fractions for the A^ 
Ac (2595)"*" ^J,~v^ and the Ag Ac (2625)"*" decays 
have been predicted, but with an uncertainty as large as 
100% g^. 

In order to reduce the systematic uncertainty on the 
final measurement coming from the branching ratios of 
these backgrounds, the following decays are searched for 
in a larger A^fiX data sample (360 pb~^): 
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A? 

Ac (2595)^ 

A+TT+TT". 



Ac (2595)+ pL-v^X where 
I]c(2455)++7r-,I]c(2455)"7r+) 



2. A° ^ Ac(2625)+Ai"i^p^ where A^ (2625)^ 

A+TT+TT". 



3. AO - 
A+^- 



Ec(2455)"7r+^-p^A: where Ec(2455)" 



4. AO Ec(2455)++7r-//-p^X where I]c(2455)+^ 
A+7r+. 



All four decays modes above contain A+7r+7r~/x~ in the 
final state. The selection criteria are the same as those 
for the A+/i^X sample (see Section llll)) . except that two 
oppositely-charged tracks are added to determine the sec- 
ondary vertex for the Ac (2595)+ and Ac (2625)+ modes, 
and one track is added for the Ec(2455) modes. In all 
cases, the pion mass is assumed for each additional track 
and each track is required to have pr > 0.4 GeV/c. The 
available four-momentum transferred to the daughters in 
the decays of these c baryons into A+ is small. Therefore, 
the mass differences Ma^^+^- — Ma^, M^^^- — Ma^, and 
MA^7r+ ~ , have a better resolution than the masses 
of the c-baryon candidates and are the figure of merit for 
detecting signal peaks. Figure [H shows the mass differ- 
ence distributions, where the numbers of signal events are 
determined by fitting the mass differences to a Gaussian 
for the signal and a kinematically-motivated line shape 
for the combinatorial background. Table IVIII summa- 
rizes the signal yields, the corresponding significances, 
and the fitted mass differences. In this table, contribu- 
tions of Ac (2595)+ in the Sc(2455)7r modes have been 
subtracted from the Ec(2455)7r modes and the signifi- 
cances of the Ec(2455)0 and Ec(2455)++ modes are com- 
bined. Systematic uncertainties on the yields are de- 
termined by varying the functions for the combinatorial 
background in the fit. This is the first observation of the 
Ac (2625)+ n^Vfj. decay. 



AO 



After estimating, with simulation, the acceptance 
times efficiency of these reconstructed decays relative to 
that of the A^ A^fi~D^ decay, and taking into account 
the false- muon background [47|, the relative branching 



ratios (i?i) are extracted: 

i3(A0^ Ac(2595)+Ai"P^ 



Ri 



A+^ i^f^ 



R2 



6 (AO 
:6 ±0. 

^(ao ^ Ac(2625)+A*-i>M 



-t-0 047 

0.126 ± 0.033 (stat) (syst) , 



B{Al 
±0 

B{Al ^ I]c(2455)07r+/i-p^ 
B{Al^ A+^i-D^) 

6(A0->Ec(2455)++^-/i-j>^) 



, +0.071 , 
0.210 ± 0.042 (stat) _^ (syst) , 



S(AO-^ Ac+/i-i>M) 

, +0.021 , 
= 0.054 ± 0.022 (stat) _^ (syst) , 

where the two Ec(2455)7r modes are averaged. Assuming 
isospin symmetry leads to the estimate: 



i?5 



6 (AO ^ Ec(2455)+7rO^-p^ 

bJa^^TT^+JFp^) 



R: 



3,4 



0.054. 



The systematic uncertainties on the relative branching 
fractions come from variation of background fitting mod- 
els and uncertainties on the low- momentum pion pT spec- 
trum and the correction to the reconstruction efficiency. 

In order to convert the above measurements of the 
relative branching fractions into absolute branching 
fractions, an estimate of B (AjJ A+/i^Pp) is required. 
A recent measurement by the DELPHI collaboration 

reported fs.O+pg (stat) +J'2(syst)^ % for this branching 

fraction [48|. However, from heavy quark symmetry, 
the semileptonic decay widths for all b hadrons are 
expected to be the same. Therefore, semileptonic 
branching fractions of the b hadrons, r''''™Yr'°*''\ 
vary only due to their lifetime differences. Since the 
AO decays to a spin-1/2 A+, contributions from both 
S and P wave amplitudes are expected. A sum of 
B{B° -> D+e-Pe) + B{B" ^ D* (2010)+ e~ De) , where 
the decays to D+ [i:»*(2010)+] correspond to the S(P) 
wave amplitudes, yield (7.33 ± 0.16) %. The number, 
7.33%, is then scaled by the world average of the ratio 
of lifetimes, tao/tbo = 0.99 + 0.10 (stat+syst) R|. The 
B (AO -> A+^-Pf,) is estimated to be (7.3 ± 0.8 ± 1.1) %, 
where the first uncertainty arises from the propagation 
of errors and the second is half of the difference be- 
tween the above estimate and the DELPHI result [49| . 
The i3 (AO A+r~i/T-) can be estimated by scaling 
S (AO ^ A+/i^i/p) by the ratio of phase space area: 
Ph.Sp. (AO ^ Atr-Dr) /Ph.Sp. (AO A^^i-P^) 
0.277. The middle portion of Table IVIIII summarizes 
the branching fractions of the Ajj semileptonic decays 
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TABLE VII: The observed number of signal events, the corresponding significance [S/y/S + -B), and the fitted mass difference 
in data for each AjJ semileptonic decay mode. 



Mode 



Yield 



31 ± 8(stat) ± 7(syst) 
53±9(stat)±5(syst) 
16 ± ll(stat) ± 7(syst) 
26 ± 12(stat) ± 9(syst) 



Significance (cr) AM [MeV/c^] 



Ac (2595)+ 

Ac (2625)+ fj.-X 

Sc(2455)°7r+/x-X 

Sc(2455)++7r-^-X 

Ec(2455) modes combined 



2.9 
5.2 



2.1 



308.47 ±0.99(stat) 
341.39 ±0.31(stat) 
166.72 ±0.69(stat) 
168.01 ±0.51(stat) 



discussed above. Uncertainties on the observed A° 
semileptonic decays and the A° Sc(2455)"'"7r°/x^P^ 
decay include uncertainties from the relative branching 
fraction measurement and uncertainty from the assumed 
K (a" ^ A+/x~i/^) . A 100% systematic uncertainty is 
also assigned to B (A° — > A+t^j/,-). 

The sum of B (AJJ — > A+/i^P^) and the branching frac- 
tions in the middle portion of Table rvTIIl is already larger 
than the inclusive A^J semileptonic branching fraction in 
the 2008 PDG summary: 

B (Ag ^ A+^-p^ anything) = 9.9 ± 2.6%. 

The following decays are, therefore, ignored in the central 
values but will be included in the systematic uncertainty: 

Ag ^ A+/o(980)/i-j>^, 

A^ — > A^7r+7r^/i^i>^(non-resonant), 

A|^ A+7r'^7r'^/i^i>^(non-rcsonant). 



TABLE VIII: The Aj semileptonic branching fractions for de- 
cays which are included in the central value (middle portion) 
and those which are not (bottom portion). All the numbers 
in parentheses are estimated uncertainties. 





Mode 


BR (%) 




A+^-j?^ 


7.3 ± (1.4) 


Ar 


-> A, (2595)+ M^i^;. 


0.9 ± (0.4) 


Ag- 


A, (2625)+ 


1.5 ± (0.6) 


Ar 


^ Ec(2455)0^+Ai-i?^ 


0.39 ±(0.23) 


Ar 


^ Ec(2455)+7rO/i-i?^ 


0.39 ±(0.23) 


Ar 


Ec(2455)++7r-/i-i>^ 


0.39 ±(0.23) 


Ag- 




2.0 ± (2.0) 


Ag- 


A+/o(980)Ai-p^ 


±(0.32) 


Ar 


A+TT+TT^ H^Vf^ 


±(0.64) 


Ar 


^ A+TT^7ry-D^ 


±(0.32) 



An estimate of these branching fractions is obtained by 
moving B (Ajj — > A^£~Di anything) upward by la. The 
remaining branching fraction is calculated to be: 



(9.9 + 2.6)% 
- 7.3% X [1 
« 1.3%. 



Modification of cr^o {pr > 6.0) /ergo (pt > 6.0) and 
estimate of B (A^ ^ A+tt") 



R1 + R2 + 3i?3 + 0.277 X B (t- 



of Ag 
Ag- 



The 1.3% is then attributed to the above decays which 
are ignored in the central value. The branching fraction 
— !■ A+7r+7r^^~i?yj is estimated to be twice that of 
A^ i:^ tt''^ ^~ V based on the isospin invariance, and 
the /o(980) mode is assumed to have the same branching 
fraction as that of the 7r°7r° mode. The bottom portion of 
Table IVIIII lists zero central values for these three decays 
and uses their estimated branching fractions above as the 
systematic uncertainties on the branching fractions. 



Equation Q requires B (A° 
ratio N. 



A+TT ) to obtain the 



G 



feed /Miad- Combining the CDF measurement 

<Tao(pt>6.0)B(A5;^A+77-) ^ , 

37l |. the world average 



'^Bo(PT>6.0)e(SO^I5 + 7r-) 

of B {B^ — > Z?+7r^), and the ratio of production cross 

^0 (pr>6.0) 



sections p \ 



PT>6.0) ' 



one may express 



A 



G 



B B^' 



(10) 
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0.28 0.30 0.32 0.34 0.36 0.38 0.40 
M(Ac7C+7C-)-M(Ac) (GeV/c') 




0.13 0.14 0.15 0.16 0.17 0.18 0.19 0.20 0.21 0.22 

M(Ac7C-)-M(Ac) (GeV/c') 




0.13 0.14 0.15 0.16 0.17 0.18 0.19 0.20 0.21 0.22 

M(Ac7C+)-M(Aj (GeV/c^) 



FIG. 9: The excited-c-baryon candidates which are as- 
sociated with a (a) Ac (2595)+ and Ac (2625) + , (b) 
Ec(2455)'', and (c) Ec(2455)++. The curves indicate fit re- 
sults to the spectra of mass differences. 



The ratio of cross sections, p, is calculated from the ex- 
pression: 

{pt > 6 GeV/c) 



(Tgo {pT > 6 GeV/c) 

CDF I 

CBuCaC 



Pt I 



(11) 



where I(^b« 
0.084 [l^. The Cbk, 



CDF I 



is the CDF I resuh: 0.236 ± 

are the correction 
factors to account for differences between the CDF I re- 
sult and this analysis in: the assumed yS(A" — > A'^i^Dg), 
kinematic acceptance, and requirements on the minimum 
Pt of AJJ and B'^. Each correction factor is explained in 
the text that follows. 

The CDF I analysis used electron-charm final states, 
such as Hb D*{2010)+e^X, Hb D+e-X, and 
Hb — > A+e^X, to measure the ratio of production cross- 
sections. The branching fraction, B (A^ A'^£~Di), was 
needed and estimated to be 7.94±0.39%, while this anal- 
ysis estimates the value to be 7.3±1.4%. The uncertainty 
1.4% is dominated by the difference from the DELPHI re- 
sult (see Section fVlI Ap . In order to be consistent within 
this analysis, a correction to the branching fraction, C^-ji, 
is applied. The value of Cs-jz is the ratio of 7.94% to 7.3% 
and found to be 1.09 ± 0.21. 

In the CDF I analysis, the A° and B'^ pT spectra mea- 
sured with fully-reconstructed decays were not available. 
In order to extract the signal acceptance and efficieiicy, 
the Nason-Dawson-Ellis (NDE) 5-quark spectrum |50 | 
followed by the Peterson fragmentation model [5l| was 
used at CDF I to obtain the pr distributions of b hadrons 
in simulation [s^. The two-track trigger allows CDF II 
to collect large samples of fully-reconstructed 6-hadron 
decays, such as A° A+tt" and B^ D+tt^, and to 
compare the pt distributions in data with those from the 
NDE-|-Peterson model. The Aj] and the 5° px spectra 
from the NDE+Pcterson model are found to be harder 
(more b hadrons at higher px) than those measured in 
data, which indicates an over-estimate of acceptance in 
the CDF I analysis, particularly for the A° decays. The 
acceptance correction factor, Cq, is the ratio of accep- 
tances using generator-level simulations with inputs from 
the measured pr spectra (identical to those described in 
Section IIV A[) and from the NDE+Peterson model: 



data— based 



/< 



,NDE+Potc 



(12) 



where an is the ratio of the kinematic acceptances of 
A|J and B'^. The value of the correction factor is found 
to be: Ca = 1-81 ^0 22 ^^r the CDF I kinematic require- 
ments [5^ . The uncertainty comes from the uncertainties 
on the measured shapes of the A" and pt distributions 
in data. 

The last correction is due to a difference in the 
minimum px requirements between the CDF I analy- 
sis [pT{Hb) > 10 GeV/c] and this analysis [pT{Hb) > 
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6 GeV/c]. By applying the same requirements to the 
A° A+7r~ and D+7r~ decays reconstructed in 

the two-track trigger data, Fig. [TOl shows that the AjJ pT 
spectrum is significantly softer (more b hadrons at lower 
Pt) than that of the j5^]. Figure fTT] illustrates the de- 
pendence of the ratio of cross-sections on the minimum 
Pt requirements, for a small [(a)] and a large [(c)] differ- 
ence between the and px spectra; the scenario in 
Fig. [TlT c) is what has been observed in data. A correc- 
tion factor, Cpj. , is required: 



Cpj, — 



_ N^o {pT > 6) / TVao (pt > 10) 

Ngo {pT > 6) / Ngo {pT > 10) 



(13) 



The Cpj, is obtained using the generator-level simulation 
with inputs from the measured px spectra of A^ and B^ 
(identical to those described in Section HV A[) . The value 
of the correction factor is found to be: Cp,, = 1.31±0.11, 
where the uncertainty also comes from the uncertainties 
on the measured pr distributions in data. 



After applying corrections C'btz^ Ca, and Cp^, p is cal- 
culated to be: 

p = 0.61 ±0.22 (CDF I) ±0.12 (DELPHI) "^H'^^ (pt), 

where the uncertainties are from the uncertainty on the 
CDF I measurement of cr^o/cr^o, the difference between 

the estimated B ^ K'^ i~ vi) for this analysis and 
that measured by DELPHI, and the uncertainties on the 
measured shapes of the A° and B^ px distributions. The 
value of p is also consistent with the result from '55']. 
The B (Ag -> A+TT") is then extracted following Eq. (fTU)) . 
with the input of the parameters listed in Table IIXI and 
is found to be: 



B{Mi ^ A+TT-) = 

0.36 ±0.07 (DELPHI) ^_^^^^^{pt) 



(syst°*''")j 



TABLE IX: Parameters for calculating B (A? A+tt") . Text 
in parentheses indicate the sources of uncertainty: the data 
sample size, general systematics, shapes of measured pr dis- 
tributions, and difference from DELPHI'S B (A° ktlTvi). 



%. 



The "syst°*'^°''" uncertainty includes the uncertainty 
on the CDF I measurement, and the uncertainty 
on G which is dominated by the world average of 
^(A+ ^pi^-7r+) [56|. This evaluation of B{kl 
A+7r~) differs from that of the Particle Data Group 
due to the differing production spectrum of the A° rel- 
ative to the B^ as described in the previous text [STf . 
The estimated value is in good agreement with the val- 
ues predicted by Leibovich et al. [13], (0.54 ± 0.18) %, 
and Cheng [ll|, (0.50 ± 0.17) %, which gives confi- 
dence in the procedure described above. The estimated 
S(A° — > A+7r~) has been used in Eq. ^ to estimate 

^othorAS 



Parameter 



Value 



G 

CDF I 

Cbtl 

Ga 

Gp^ 

B(B" D+TT- 



0.82 ± 0.25(stat ® syst) ± 0.06(pt) 
0.236 ±0.084(statesyst) 
1.09±0.21(DELPHI) 



1.81 



-0.22 



(pt) 



1.31±0.11(pt) 
(0.268 ±0.013)% 
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Exp. const = 2.41 ±0.61 

Exp. slope = 0.077 ± 0.020 
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FIG. 10: Comparison of the reconstructed A^ and 5° pr 
spectra in data. The negative slope (3-4(t away from zero) of 
the ratio of A^J to B" histograms indicates that the pr (A°) 
distribution is softer (more b hadrons at lower pr) than the 
Pt {B^^) distribution. In order to have a fair comparison of 
Pt spectra, the same requirements are applied to the B'^ and 
Aj candidates [Hj, while Fig.[3]has different selections for the 
B" and A^. Nevertheless, the pt spectra, used as inputs for 
the correction factors Ca and Cpj,, have been corrected for 
acceptance and efficiency and are identical to those described 
in Section ITV^ 



feed 



V^had (see Section lyiB]) . 
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FIG. 11: Examples of the AjJ and -B" pr spectra [(a), (c)] and the dependence of the production cross-section ratio on the 
minimum pt requirements, p^™ [(b), (d)]. Figures [TH a') and (b) show the case where both hadrons have similar spectra; the 
ratio of the integrated areas underneath the spectra, from p^™ and above, depends little on the value of p^™ . Figures flTT c) 
and (d) show that the Aj pt spectrum is significantly softer (more h hadrons at lower pt) than the pr spectrum; the ratio 
of the integrated areas depends strongly on the value of p^^^ . 



VIII. SYSTEMATIC UNCERTAINTIES 

The relative branching fractions {R), with statisti- 
cal uncertainty only, can now be extracted from Eq. ([T]): 



R = 



g(Ag^A+/i-z7,) 



= 16.6 ± 3.0 (stat) . 



A check of internal consistency is performed by dividing 
the data and simulation samples into several groups of 
independent subsets, according to the time period, ver- 
tex position, ct and pr of the A+ candidate, ct* and pr 
of the A° candidate, etc. Figure fT2l shows that the R of 
each subset for each group is consistent with those of the 
other subsets in the same group. The result of this check 
also proves that there is no major problem in the detec- 
tor, trigger, reconstruction, or simulation which produces 
bias in the measurement. 

The systematic uncertainties on R may be classified 
as internal and external. Internal uncertainties are those 



which affect the final measurement through their effects 
on the observed yields, the numbers of false-muon and 
bb/cc background events, and the modeling of acceptance 
times efficiency. External uncertainties are those from 
production fractions and branching ratios which are used 
in Eqs. ^ and (|Aip to determine iVfeod- The input value 
for each systematic source is varied by ±lcr where a is the 
uncertainty on the input value. The resulting difference 
in R from the central value is the systematic uncertainty. 
The following text describes how the uncertainty for each 
systematic source is obtained. 



A. Internal systematic uncertainties 

The signal yields, A^had and A^incisemi, are affected by 
the background functions which describe the mass spec- 
tra of misidentified or partially-reconstructed decays of 
b hadrons. The systematic uncertainty on the M^^^ fit- 
ting model is estimated by changing the relative fraction 
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B(A°^AJH-V^)/B(A°^A,-7i-) 



o o 


\3 GJ 
D O 


Time period 1 
Time period 2 
Time period 3 
Time period 4 


J.. . . 1 . 


1 1 1 1 1 1 1 1 


\ks vertex z > cm 
\ks vertex z < cm 






PT(Vs)>-'2GeV/c 
PT(Vs)<''2GeV/c 








Pt(A^) > 8 GeV/c 
Pt(A*) < 8 GeV/c 






ct*(A°)>6.64cm 
ct*(A^) < 0.04 cm 




ct(A*) > 0.005 cm 
ct(A;) < 0.005 cm 






2<x^(4.J<15 

Xr^,{\J < 2 






6.5<X?,(AJ)<14 
X?j(Aj)<0.5 








4.3 < IW(A^n) < 5.7 GeV/c^ 
■3.7 < IVl(A<,n) < 4.3 GeV/c^ 









FIG. 12: Internal-consistency check of the A° relative branch- 
ing fractions. The uncertainty on each point is statistical only. 
Each independent group is separated by a horizontal dashed 
line. The solid bands indicate the relative branching frac- 
tions with their statistical uncertainties from the complete, 
undivided samples. 



of the contributing decays in each background function. 
The fragmentation fraction and the branching ratio of ev- 
ery contributing decay are varied independently accord- 
ing to their uncertainties by ±la [58]. After combin- 
ing these contributing decays according to their modified 
fragmentation fractions and branching ratios and pro- 
ducing a new simulated mass spectrum, the parameter 
values for each background function are re-determined. 
The Ma^tt spectrum in data is re-fit again using the new 
background function and the systematic uncertainty is 
taken from the deviation of the A+7r~ yield from 

the central value. Correlations between different param- 
eters have been taken into account. For the fitting of 
MpKTn since no branching ratio assumptions are made, 
no systematic uncertainty is assigned. 

The uncertainty on the false-muon estimate is driven 
by: 1) the size of the sample used to measure the false- 
muon probability, 2) the fit to the weighted Ac mass dis- 
tributions, and 3) the probability of a hadron track being 
a p, a K, or a tt, which is determined from simulation. 
The resulting changes in the number of false-muon events 
from the three sources above are added in quadrature 
and already listed in Table IIVI The size of the bb/cc 



background contribution has a 100% systematic uncer- 
tainty, due to a lack of knowledge of the relative bb/cc 
production rates between different processes [s^ and the 
10-50% discrepancy of the inclusive hadron production 
cross-section between PYTHIA and data (Appendix [C| . 

The uncertainty on the modeling of acceptance times 
efficiency for signal and background processes arises 
from: the size of simulation samples, the shapes of the 
measured A° and B° pt spectra, the efficiency scale fac- 
tors of the AfJ semileptonic decay model/muon recon- 
struction/XFT trigger (see Section HV|) . the amount of 
material in the detector simulation, the A^ lifetime, the 
A+ Dalitz structure, and the and A+ polarizations. 
The A° and B° px distributions used as inputs for BGEN- 
ERATOR are varied according to the uncertainties on the 
exponential slopes of data to MC ratios (shown in Fig. [3]). 
The uncertainties on the efficiency scale factors for the 
A^ semileptonic decay model/muon reconstruction/XFT 
trigger are listed in Table [TTl The uncertainty from 
the detector material is obtained by switching off the 
hadronic interaction in the detector simulation and mul- 
tiplying the efficiency difference from the central value by 
25%. The 25% is a quadrature sum of the 15% underes- 
timate in the amount of material and the 20% difference 
between the GHEISHA and FLUKA models [H, 113 . 
The A° lifetime used as an input for bgenerator is 
varied according to the uncertainty on the world aver- 
age 7] . The effect of the A+ Dalitz structure is studied 
by varying branching fractions of the resonant and non- 
resonant A+ pK'TT'^ decays measured by E791 ^ by 
their uncertainties. The unpolarized A° and A+ simula- 
tion samples have been used to obtain the central values 
of acceptance times efficiency of A° decays. For the sys- 
tematics study, angular distributions in simulation are re- 
weighted according to all combinations of the A° and A+ 
polarization states: ±1, assuming the extreme scenario 
where the A" and A+ baryons are 100% polarized. The 
difference in the kinematic acceptances between the sim- 
ulation with re-weighted angular distributions and the 
simulation with unpolarized Aj| and A^ is used to assign 
a systematic uncertainty on R. 



B. External systematic uncertainties 

There are two types of external systematic uncertain- 
ties. The first type is denoted as the "PDG" uncer- 
tainty and includes uncertainties on: the world average of 
B {^B^ £)+7r~), the CDF I measurement of aj^o/ago, 

<y^o (pT>6.0)e(A5;-*A+7r-) 

the CDF measurement of — —, — ^„ nMg/an — ?t+ — and 
the measured branching fractions of the four new A^ 
semileptonic decays relative to that of the A" A+/x~P^ 
decay (see Section IVII A[) . The second type is denoted 
as the "EBR" uncertainty and comes from unmeasured 
branching fractions estimated from theory. A 5% uncer- 
tainty is assigned to the estimated branching fractions 
of the excited c-hadron decays [i^. A 100% uncertainty 
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is assigned to the other unobserved 6-hadron decays to 
cover the wide range of theoretical predictions [HJ . Note 
that the uncertainty on the estimated B (A" A+^~P^) 
does not affect the final measurement because it affects 
the branching fractions of A° — > A+tt^ and other A° 
semileptonic decays in the same way so that any change 
completely cancels. 

C. Summary 

The fractional systematic and statistical uncertain- 
ties on the A^ relative branching fraction are summa- 
rized in Table |Xl The leading sources of internal sys- 
tematic uncertainty are the mass fitting model, the 
shapes of the measured pr spectra, and the A^ semilep- 
tonic decay model. The PDG uncertainty is domi- 
nated by the world average of ;B (A+ ^ pK~Tr'^) and 
the CDF I measurement of cr^o /cr^n which have been 

used to extract B (A[J A+tt^) 62]. The uncertainty on 
B (A+ — > pK^TT^) may be reduced in the near future by 
more precise measurements proposed by Dunietz [63j and 
Migliozzi [g^] . The EBR uncertainty is dominated by the 
branching fractions of non-resonant A^ —> K'^ t:'^ it~ ^~ v ^ 
and A+TT^Tr'^/i^t'^ decays. The A" relative branch- 

ing fraction with complete uncertainties is found to be: 

g(A°^A+M-z>^) _ 

16.6 ± 3.0 (stat) ± 1.0 (syst) "^^'^(PDC) ± 0.3(EBR). 

The uncertainties are from statistics (stat), internal sys- 
tematics (syst), world averages of measurements pub- 
lished by the Particle Data Group or subsidiary measure- 
ments in this analysis (PDG), and unmeasured branching 
fractions estimated from theory (EBR), respectively. 

IX. MEASUREMENTS OF THE B° RELATIVE 
BRANCHING FRACTIONS 

The same analysis technique used for the A^ samples is 
applied to the 5° decays. This section only describes the 
difference in the details of event reconstruction, yield de- 
termination, and background estimation and summarizes 
the systematic uncertainties. 

A. Reconstruction of the B'^ candidates 

The following decay modes are reconstructed in the 
data collected with the two-track trigger: 

1. D+TT- and D+^-X, where D+ 

2. 5° ^ D*(2010)+7r- and (2010)+^"^= where 
D*(2010)+ ^ D°t:+,D^ K-tt+. 



TABLE X: Summary of statistical and systematic uncertain- 
ties for the A^ mode. The an is the uncertainty on the A|J 
relative branching fraction, R. 



Source 




Mass fitting 


-1-3.8 
-3.1 


False /i 


1.0 


bb/cc background 


0.3 


Simulation sample size 


2.0 


6-hadron px spectrum 


-1-0.0 
-2.9 


A° decay model 


3.3 


XFT/CMU efficiency scale factor 


0.4 


detector material 


1.3 


A^ lifetime 


0.3 


A+ Dalitz 


0.4 


A|J, A+ polarizations 


2.2 


Sum of internal 


6.3 


PDG 


+ 15.6 
-20.4 


Estimated branching fractions 


2.1 


Statistical 


17.8 



The requirements on the and A° candidates are kept 
as similar as possible. 

For the reconstruction of K^n^n^ decays, the 

pion mass is assigned to the two positively charged tracks 
and the kaon mass to the negatively charged track. The 
invariant mass of the three tracks (Mif^jr), as computed 
by a three-track kinematic fit, is required to be in the 
range 1.74 - 2.00 GeV/c^. The D*{2QIQ)+ signals are 
reconstructed by first looking for D'^ K^tt^ candi- 
dates. A two-track kinematic fit determines the D'^ ver- 
tex, and the invariant mass of the two tracks (Mktt) is 
required to be within the range 1.820 - 1.906 GeV/c^. 
Then, the pion mass is assigned to an additional pos- 
itively charged track. This third track is expected to 
have a low px due to the small four-momentum transfer 
in the D*(2010)+ — > Z3°7r+ decay. However, a minimum 
Pt requirement of 0.5 GeV/c is imposed to ensure a good 
measurement of the pion track. For the £)*(2010)+ can- 
didate, the mass difference, MKTrTr—MK-n, must be within 
the range 0.14 - 0.18 GeV/c^. 

In order to form a candidate, the £>+ and 
D*(2010)+ candidates are then combined with an ad- 
ditional negatively charged track which satisfies the re- 
quirements described in Section IIIII After the four- 
track kinematic fit, the values of Mk-k-r for the D+ and 
Mktttt — Mktt for the Z)*(2010)"'" must be in the range: 
1.8517 - 1.8837 GeV/c^ and 0.143 - 0.148 GeV/c^ for 
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the hadronic candidates; 1.74 - 2.00 GeV/c^ and 0.14 - 
0.18 GeV/c^ for the inclusive semileptonic candidates. 
The four-track invariant mass, Mjj^ and Md*^, must 
be within 3.0 - 5.3 GeV/c^ for the semileptonic decays. 
Selection criteria for the following variables: pt of the 
fourth iJ^-candidate track [pr (tt^, ^^)], px of pr of 
D*(2010)+, and combined pr of the four-track system, 
of the and Z?" vertex fits, and of the four- 
track kinematic fits, ct of the and D'^ candidates, 
and ct* of the candidate, are also optimized using the 
simulation and data of hadronic modes, as described for 
the A° sample. Table IXll lists the optimized values. 



B. Determination of the B° yields 

Figure [TSf a) shows the fit result for the Mj:,^ spec- 
trum. The 5° D^TT^ yield returned by the fit is 
579 ±30. The signal peak at Md^ « 5.3 GeV/c^ and the 
combinatorial background are described by a Gaussian 
function and an exponential, respectively. The ratio of 
the number of doubly Cabibbo-suppressed decays rela- 
tive to that of the signal mode, Ndk/Nd-k, is Gaussian- 
constrained to the value for the relative branching ratio 
from the PDG, convoluted with the efficiency from the 
full simulation. The constrained value is 0.073 ± 0.023. 
Backgrounds from the other 6-hadron decays consist of 
the following decays. The Dfn' decays, where 

D+ (/)(1020)7r+, 0(1020) -> K+K~ and the pion mass 
is assigned to one of the kaons, appear as a peak at 
around 5.31 GeV/c^. Misreconstructed A+tt" de- 

cays, where A+ pK~TT~^ and the pion mass is assigned 
to the proton, form a broad peak around 5.4 GeV/c^. 
The backgrounds from the B° D+K', B° D+tt", 
and A(] —> A+tt" decays are combined and indicated 
by the black filled region. The i?° D^p~ decays, 
where p~ — > tt'^tt" and the tt*' is not reconstructed in 
the event, have a triangular mass distribution which 
peaks at « 5.1 GeV/c^. The B° D*{2010)+n- de- 
cays, where _D*(2010)+ D+tt'^ and the 7r° is not re- 
constructed, have a double-peak structure. This struc- 
ture is consistent with the spin-1 £'*(2010)+ being po- 
larized. This polarization results in the 7r° from the 
£'*(2010)"'" decay having a momentum preferentially par- 
allel or anti-parallel to the momentum of £'*(2010)+. 
The B° D+p~ and B^ D*{2Q\Q)+t:- backgrounds 
are combined and indicated by the dark-gray filled re- 
gion. The remaining partially-reconstructed decays of 
h hadrons, H\j — > D+X, have a monotonically falling 
distribution (hatched region). The determination of the 
background shapes and the estimation of systematic un- 
certainty are similar to those in the A^ system. 

Figure [T^Tb) shows the fit result for the Mk-k-k spec- 
trum for events with muons. The inclusive D'^ p^ X yield 
returned by the fit is 4720 ± 100. The signal peak at 
Mk-ktt ~ 1.9 GeV/c^ is described by a Gaussian function. 
The combinatorial background (light-gray filled region) is 
parameterized by a first-order polynomial. Misidentified 



Hi, Dfp~X decays (black filled region), where the 
mass of at least one Df daughter has been misassigned, 
appear in the mass window of interest. The dominant 
contributing D+ decay modes are 0(lO2O)7r"'", 
D+ K*°K+, and D+ ^non-resonant K+K-t:+. The 
function parameters describing the shape of misidenti- 
fied Df spectrum are obtained from the B'^ — > Df p,~v^ 
simulation; in this simulation, the D'^^s are forced to de- 
cay only to the final states which can yield misidentified 
mass in the Mx-mT window. The number of the D'^ back- 
ground events is constrained to the estimated number of 
Hb Dfp,~X events in the data (iV^^^) as described 
below. First, the following decay mode is reconstructed 
in the data: Hb Dfp~X, where Df -> 0(lO2O)7r+ 
and 0(1020) K^K~ . The narrow 0(1020) resonance 
provides a good handle for removing the combinato- 
rial background of D+. Second, the fraction from the 
0(lO2O)7r"'" mode relative to all contributing £)+ decays 
(i?0^) is extracted using the world averages of the 
branching ratios, and the acceptance times efficiency de- 
termined from the full simulation. Then, Nd^^ is simply 
the yield of the 0(lO2O)7r+/Li~ mode in data divided by 
i?07r- The value of No^fi for the constraint is 1812 ± 160. 
The systematic uncertainty is assigned by independently 
varying the ratio of the branching fraction of one specific 
decay relative to that of the £)+ 0(1020)77+ decay 
by zblfT, since the branching fractions of all Dg decays 
have been measured relative to that of the 0(lO2O)7r+ 
mode Q. 

Figure [TSTc) shows the fit result for the Md'tt spec- 
trum. The B° D*(2010)+7r- yield returned by the fit 
is 106 ± 11. The analysis of the D*(2010)+7r- signal and 
backgrounds is similar to that in the D~^Tr~ mode. The 
only difference is that extra constraints are imposed due 
to the small size of the D*(2010)"'"7r~ sample. The width 
of the signal Gaussian (Jd'tti the ratio , and the ra- 

tio of backgrounds No'p/ {No'p + NH.^rcmaining d-x), 

are constrained to 0.0259±0.0012 GeV/c^, 0.071±0.019, 
and 0.242 ± 0.008, respectively. The systematic uncer- 
tainty is assessed in the same way as in the A|J and the 
ZJ+TT" modes. 

Figure [T^Td) shows the fit result for the Mktttt — Mk-k 
spectrum for events with muons. The likelihood fit 
for the D*(2010)"'"/i~Ar mode is performed in the mass 
window Mk-k-k - Mktt = 0.14 - 0.18 GeV/c^ whereas 
Fig. 113r d) shows a more restricted mass range near the 
signal peaks. The inclusive D*{2{)\Q)^p~X yield re- 
turned by the fit is 1059 ± 33. The signal peak at 
Mk-ktt — Mktt ~ 0.145 GeV/c^ is modeled by two Gaus- 
sian distributions with a common mean and different 
widths. The combinatorial background (light-gray filled 
region) is parameterized by a constant, while the back- 
ground from other 6-hadron decays with misidentified c- 
hadron daughters is found to be negligible. The size of 
the combinatorial background is very small due to the 
requirement that is consistent with the world av- 

erage D'^ mass, the minimum requirement on the mass 
Mo'fi, and the minimum requirements on the px and the 
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TABLE XI: Optimized requirements for reconstructing the B° D+tt', D+fi'X, B° D*(2010)+7r~, and D*(2010) + /x"X 
decays. 





■> D+Ti- 


5° L>*(2010)+7r- 
D*(2010)+Ai-X 


Vt (tt^ , u,~) 

1-' 1 \" 1 J 


> 2 GeV/c 


PT(7r ,^1 ) 


> 2 GeV/c 


PT {D+) 


> 5 GeV/c 


Pt(^*(2010)+) 


> 5 GeV/c 


PT (4trks) 


> 6 GeV/c 


PT (4trks) 


> 6 GeV/c 




< 14 


Xl, {D°) 


< 16 


Xl^, (4trks) 


< 15 


Xl^, (4trks) 


< 17 


ct{D+) 


> —30 /im 




> -70 /im 


ct* {&) 


> 200 Aim 


ct* {B°) 


> 200 ^m 



number of SVX hits for the low-momentum pion from the 
£)*(2010)+ decay (Section |TX2|. The fitting function for 
this spectrum does not use any branching ratios and no 
systematic uncertainty is assigned. 

Table IXIII summarizes the i?" hadronic and inclusive 
semileptonic yields and the probability of correspond- 
ing fits. Each model describes the data well, as indicated 
by the probability. 



C. Compositions of tlie inclusive semileptonic data 

The procedures for estimating the backgrounds to the 
B" -> and B" -> D* (2010)+ fi'Df, decays are 

similar to those described in Section IVII The follow- 
ing describes the differences when estimating the feed-in 
background, iVfced, in the B° system. Unlike the A° sys- 
tem, many decays of b and c mesons have been measured 
by other experiments [13, HH, [H, HI] , and serve as inputs 
to the EVTGEN decay package. In addition, evtgen also 
includes estimates of branching fractions for decay modes 
which have not yet been measured. Therefore, all possi- 
ble decays which may contribute to the A^food in the 
control samples are studied using the PDG summary and 
the default evtgen decay table 0, HI . 

The feed-in backgrounds to the B'^ fi^v^^ and 

B^ L)*(2010)"'"/i~P^ signals fall into two categories: 

1. Semileptonic decays of B^/B^ /Bl, which in- 
clude either additional particles {e.g., B^ — > 
D^n'^fi^Vfj^) or a higher mass c meson with sub- 
sequent decay into the c-meson signal {e.g., B^ — > 
i:i*(2010)V"f^M' £'*(2010)+ D+ttO) 

2. Hadronic decays of b mesons into two c mesons: one 
c meson decays hadronically in a reconstructed fi- 
nal state, the other c meson decays semileptonically 
{e.g., B° -> D+DJ, Dj ^ 0(1020)^"^^^). 



Branching fractions of the B DD decays relative to 
the signal are all below 3%. A generator-level study indi- 
cates that they are further suppressed after a minimum 
requirement on the four-track invariant mass, M/j^^j*)^, 
and therefore, contribute less than 1% to the signal. 
Backgrounds from b mesons decaying semileptonically to 
more particles or higher mass c mesons are also reduced 
or eliminated by the same minimum mass requirement. 

Tables IXIIII IXTVl summarize the feed-in backgrounds 
which contribute > 1% to the i?° ^~v^ and the 

— > L'*(2010)+^~P^ decays. The definition of quanti- 
ties listed in each column follows Table |Vl Only these 
decays are subtracted from the inclusive semileptonic 
yield. The leading background to B^ D^^~D^ is 
ijo ^ D*(2010)+Ai-i>^ where D*{2QIQ)+ D+tt°. The 
leading background to B° D*{2010)+n-D^, is B- 
L»i(2420)V~i^M where £)i(2420)" ^ i:»*(2010)+7r-. 
Combining information compiled in the PDG, back- 
grounds from B D'^*'**^ fi^ Df^X which are not consid- 
ered in Tables IxnTHxTVl contribute less than 2% to the 
signal. The estimates of A^faise/i and Nf^j, for the i?" are 
identical to those for the A°. Table IXVl lists the results. 
Figure [HI shows the Mktttt and Mk-kit — Mk-k distribu- 
tions weighted with muon-misidentification probabilities 
and the results of the fit. 

The compositions of the inclusive ^~ and 
£'*(2010)+/L(~ samples are summarized in Table IXVII 
The dominant signal contamination is from the feed-in 
background. The second largest background arises from 
false muons. The smallest background source is from 
66/cc. 

D. Systematic uncertainties 

Figure \TE\ shows a summary of internal consis- 
tency checks. The fractional systematic and statis- 
tical uncertainties on the relative branching frac- 
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TABLE XII: Observed number of events in each decay mode determined from the unbinned, extended likelihood fit, x^/NDF, 
and the corresponding probability computed to indicate quality of the fit. 



Mode Yield xVNDF Prob (%) 

B" ^ D+TT- 579 ± 30 80/91 78.9 

D+^i-X 4720 ± 100 47/31 3.40 

i?o ^ i:)*(2010)+7r- 106 ±11 21/12 5.40 

D* (2010)+ fi- X 1059 ±33 108/93 14.1 
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FIG. 13: Results (curve) of the unbinned, extended likelihood fits for determining the numbers of 5° candidates: the hadronic 
modes (a) Ado-n and (c) Md'tt, and inclusive semileptonic modes (b) Mktvtv and (d) Mk-rtt ~ Mk-w- The filled histograms 
indicate various backgrounds. 
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TABLE XIII: The feed-in backgrounds to the D+^i~u^ signal. For the ^ D+K°fi~Uf, decay, the 

T{B°)/ [r{B~) +r{B°)] from the PDG is used to obtain iVfocd/A^had [H, E^. Numbers in parentheses are estimated un- 
certainties for the unmeasured branching fractions [63] • The definition of quantities listed in each column follows Table FVl 



Mode 


B (%) 






( Nto^i \ 


■i 

event 




V A'had ) 




^ D+n- 


0.268 ±0.013 


1.000 






579 ± 30 












4720 ±100 




2.17±0.12 


0.455 ±0.004 


3.688 


1.000 





fiO D*(2010)+Ai-P^ 5.16±0.11 0.372 ±0.004 2.314 0.627 1340± 69 

^D+TT°/j 32.30±0.64 

B" D+n^fi-Pf, 0.30±(0.30) 0.165 ± 0.003 0.185 0.050 107± 6 

B"^D+T~P^ 1.00±0.40 0.100±0.004 0.065 0.018 37± 2 



^ ii'v^.vr 17.36 ±0.05 

B- ^ Di{2mf^i-D^ 0.40±0.07 0.278±0.005 0.134 0.036 78± 4 
^ i:>*(2010)+7r- 

-^L>+7r"/7 32.30±0.64 

B- ^ D[{2AiQfiJL-Vf, 0.37±(0.37) 0.273 ±0.005 0.081 0.022 47± 2 

L>*(2010)+7r- 66.67±(3.33) 

-^L>+7r"/7 32.30±0.64 

B- ^ D+TT-fi-Df, 0.42±0.05 0.165±0.003 0.259 0.070 150± 8 

B° D+K^fi-Di, 0.30 ±(0.30) 0.137 ±0.004 0.064 0.017 37± 2 



TABLE XIV: The feed-in backgrounds to the B° D* (2010)+ n~ signal. 



Mode B (%) ; '-^ (^^^^) N,. 





^ i:>*(2010)+7r- 


0.276±0.013 


1.000 






106± 11 


D* (2010)+/!" a: 










1059 ±33 


& - 


i:>*(2010)+Ai"i^^ 


5.16±0.11 


0.447 ±0.006 


8.361 


1.000 




B°- 




0.81 ±(0.32) 


0.349 ±0.008 


0.341 


0.041 


36± 4 




^ D*(2010)+7r0 


33.33±(1.67) 










B°- 




0.37±(0.37) 


0.336 ±0.008 


0.150 


0.018 


16± 2 




^ D*(2010)+7rO 


33.33±(1.67) 












D*(2010)+7rV"t^M 


0.10±(0.10) 


0.239 ±0.006 


0.086 


0.010 


9± 1 


B"- 


D*(2010)+T-j>^ 


1.60±0.50 


0.136 ±0.005 


0.137 


0.016 


14± 2 






17.36±0.05 










B- 


Di{2420)y-D^ 


0.40 ±0.07 


0.356 ±0.008 


0.516 


0.062 


55 ± 6 




^ D*(2010)+7r- 












B- - 


D[{2430ffi-P^ 


0.37±(0.37) 


0.351 ±0.008 


0.314 


0.038 


33± 3 




^ L>*(2010)+7r- 


66.67±(3.33) 










B- 


D*{2Q1Q)+Tr- fj.-Df, 


0.61 ±0.06 


0.242 ±0.006 


0.534 


0.064 


57± 6 



tions are summarized in Table IXVIII The leading 
sources of internal systematic uncertainties are the 
mass fitting models and the shape of the measured 
B^ Pt spectrum. The dominant PDG uncertain- 
ties come from B (5° D+7r~) for the Z?+ mode, 
and B (5° i:»i(2420)+^-p^) for the D*(2010)+ mode. 
The dominant uncertainties on the estimated branching 
fractions come from B {B^ D'^ ■k^ fiT T> ^ for the D"*" 



mode and B {B- D[{2430)"pi~Pf,) for the D*(2010)+ 
mode. The B'^ relative branching fractions with complete 
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FIG. 14: The invariant mass distributions produced from data with a hadron track and a c- meson candidate in the 

final state, after weighting the hadron track with an average muon-misidentification probability (Pavg): (a) Mktvtv and (b) 
Mktvtt — Mktt- The curves indicate the results of the fil- 



TABLE XV: The estimated sizes of the false-muon, bh, cc, 
and feed-in background contribution to the — > ^~v^^ 
and B° -* _D*(2010)+ ff, signals and the observed yields 
in data. Uncertainties on the bh, cc, and feed-in backgrounds 
are statistical only. 







D*{2010)+n 


^falsc^ 


230 ±19 


44 ±3 




0.08±0.01 


0.08±0.01 


A^cc/iVhad 


0.05±0.01 


0.05±0.01 


A^fecd/A^had 


3.10±0.03 


2.08±0.02 


A^had 


579 ±30 


106±11 


inclscmi 


4720 ±100 


1059 ±33 



TABLE XVI: The composition of the inclusive D'^ fi and 
D*(2010) + ^i" data sample. 



semi (%) 

Signal 55.5 73.7 

False muon 4.9 4.2 

bb/cc 1.6 1.3 

Fccd-in 38.0 20.8 



ments in this analysis (PDG), and unmeasured branching 
fractions estimated from theory (EBR), respectively. 



uncertainties are found to be: 

9.9 ± 1.0 (stat) ± 0.6 (syst) ± 0.4(PDG) ± 0.5(EBR), 
^ £)*(2010)+M-i?p) _ 
6(50 i:>*(2010)+7r-) ~ 
16.5 ± 2.3 (stat) ± 0.6 (syst) ± 0.5(PDG) ± 0.8(EBR). 

The uncertainties are from statistics (stat), internal sys- 
tematics (syst), world averages of measurements pub- 
lished by the Particle Data Group or subsidiary measure- 



X. RESULTS 

The and 13^ relative branching fractions are mea- 
sured to be: 

B (Ag ^ A+m-^m) _ 
B{Kl^Kt^-) 

16.6 ± 3.0 (stat) ± 1.0 (syst) ^^'^(PDG) ± 0.3(EBR), 
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TABLE XVII: Summary of statistical and systematic uncertainties for the B modes. The an is uncertainty on the B relative 
branching fraction, R. 



Source 




if (%) 

Z3(s''^Z5*(2010) + m"Pm) 


rif DO ^7-)-|-__\ 

O I -D - — > U'TT 1 


\ 


^Ll*(2010)+7r-) 


Mass fitting 


4.1 




< 0.1 


Fcilse jj, 


0.7 




0.4 




2.9 




1.7 


Rimnlrition sj^mnlp si^p 

k-/i.xxx i^xcaj u xwxx o capxxx lj±v_> oxzj 


1.6 




1.7 


px spectrum. 


3.0 




2.4 


jLb i /CMU emciency scale lactor 


0.5 




0.4 


detector material 


1.7 




1.3 


Sum of internal 


6.3 




3.7 


PDG 


4.1 




2.8 


Estimated branching fractions 


4.7 




4.9 


Statistical 


9.7 




14.1 



B ^ D+TT-) 
9.9 ± 1.0 (stat) ± 0.6 (syst) ± 0.4(PDG) ± 0.5(EBR), 

B{B° ^ D*(2010)+7r-) ^ 
16.5 ± 2.3 (stat) ± 0.6 (syst) ± 0.5(PDG) ± 0.8(EBR). 

The uncertainties are from statistics (stat), internal sys- 
tematics (syst), world averages of measurements pub- 
lished by the Particle Data Group or subsidiary measure- 
ments in this analysis (PDG), and unmeasured branch- 
ing fractions estimated from theory (EBR), respectively. 
The control sample results are consistent with the ra- 
tios published by the 2008 P DG ?| at the 1.3cr and 0.8cr 
level, respectively (see Table IXVlIIj) . The measured ra- 
tio of AJJ branching fractions is compared with the pre- 
dicted value based on HQET. The prediction has a « 30% 
uncertainty and is obtained by combinin g th e results of 
Huang et al. [H and Leibovich et al. p^|68|. Figure [TCI 
shows the consistency between this measurement and the 
theoretical prediction. 

XI. CONCLUSION 

Using data from an integrated luminosity of 
« 172 pb"^ collected with the CDF II detector, 1237 ±97 
A+fi~X and 179 ± 19 Aj^ — > A+tt" signal events are 
The branching fractions of the four new A° semilep- reconstructed. The large A" sample enables the measure- 
tonic decays relative to that of A° A+/i~i/^ are mea- ment of B (A" A+/i~i/^) /B (A|J A+7r~) and the 



sured to be: 

b(^A°^ A,{2595)+ fi-Df, 



+0 047 

0.126 ± 0.033 (stat) _^ (syst) , 



.126 d 

6(aO-> Ae(2625)+M"P^ 
S(AO^A+^-p^) 

0.210 ± 0.042 (stat) (syst) , 

^ ' -0.050 ^ ' 



1 r g(Ag^S,(2455)"7rV-p^) 
2[ B{Al^A+^~v^) 

6(Ag^S,(2455)++^-Ai-i^^) 



B{Al^A+^i-v,) 

0.054 ± 0.022 (stat) ^^f'^^^ (syst) 
— 0.018 
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TABLE XVIII: The B° relative branching fractions measured in this analysis and those published in the 2008 PDG 0]- The 
measurements of this analysis include both the statistical and the systematic uncertainties. 



Mode 



B{B°^D' (2010) + fi-u^) 
B(_B»^D'(2010) + 7r-) 



PDG 



This Analysis 



8.1 ±0.6 9.9 ±1.3 
18.7±1.0 16.5±2.6 



comparison to the predictions of heavy quark effective 
theory. The uncertainty is dominated by the size of the 
data sample, the world average of B (A+ pK~TT^), and 



the CDF I measurement of cr^o / ago ■ 
trol modes, B (5° ^ D+^-D^) /B 



Ratios for the con- 
^ D+TT-) and 



B (B° ^ D*{2010)+fi-Df,) /B {B° i:>*(2010)- 



are found to be in good agreement with the world 
averages [7]. For the first time, the semilep- 
tonic decay A° 
observed and three 
Ac(2595)+/i"i 



Ac (2625) /i z/^ has been 



other 



semileptonic decays. 



and AJ; Ec(2455)"''^7r~/i~P^, have been recon- 

structed, using data from an integrated luminosity of 
« 360 pb~^. Measurements of the ratios of their branch- 
ing fractions to the branching fraction of A° A+/x~Pp 
have been performed. Finally, the transverse-momentum 
distribution of the Ag baryon produced in pp collisions 
is found to be softer (more b hadrons at lower pt) than 
that of the meson; this results in a new estimate for 
B (A" A+7r~) in better agreement with the theory 
than the PDG evaluation. 
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APPENDIX A: SEMILEPTONIC DECAYS OF b 
MESONS TO BARYONS 



The number of feed-in events from semileptonic decays 
of b mesons to baryons {Ni^°l°'^) is also normalized to the 
Ag A+TT^ yield in data (A'had) and has an expression 
similar to that of Eq. ^ : 



Armcson 

t^-B„ aAPT> 6.0) 



ctao(pt>6.0) B{Al^A+TT-)e 



-(Al) 



where ctb^ ^ {pT > 6.0) and a^^o {pT > 6.0) are the pro- 
duction cross-sections of b mesons and Aj] baryons for px 
greater than 6 GeV/c. 

A list of 6-meson decays that may contribute to the 
A^/i~A" sample is obtained from an inclusive sample 
of 6- meson semi-muonic decays generated using pythia. 
After applying the trigger and analysis requirements to 
the PYTHIA generated events, the maximum contribut- 
ing decays are found to be A+n/^~i>^ and B~ 



While there are measurements of branching 



ratios of the 6-meson hadronic decays to baryons, e.g., 
B^ — > Ajj"p7r~'"7r~, there is only an upper limit for the 
semileptonic decay of B^: 

B {B- A+pe-j>e) < 0.15%. 

Assuming the branching fractions of the muon-neutron 
and muon-proton final states are the same as that of the 
proton-electron final state, the value of this upper limit 
is then taken for the branching fraction of the B'^ 
A+n/i^j/p and the B^ — > A^pfi~v^ decays. The ratio 



N 



for example, is then given by: 



N 



AO-^A+TT- 



aso {pT > 6.0) 13 {B" A+n^-p^.) ego^A+ftM-p^ 
^Ao {pT > 6.0) B (Ag ^ A+TT-) eAO_.A+.~ 

1 0.15% ega^A+ 



G B B^^D+TT^ 



(A2) 



>A+7r- 
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FIG. 15: Internal-consistency checks of the relative branching 
fractions measured in the two control samples: (a) — ^ 
D+ decays and (b) B° D*(2010)+ decays. The uncertainty 
on each point is statistical only. Each independent group 
is separated by a horizontal dashed line. The solid bands 
indicate the relative branching fractions with their statistical 
uncertainties from the complete, undivided samples. 
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FIG. 16: Comparison of the measured 

B {A° A+fj.'u^) /B (Aj A+tt") with the theoreti- 
cal prediction based on HQET by Huang et al. Q and 
Leibovich et al. [T^]. The measurement includes both the 
statistical and the systematic uncertainties. 



where G is the CDF measurement [37| : 

aAo(pT>6.0)6(Ag->A+^-) 



G = 



cTsoipT > 6.0)B{B° ^ D+TT- 



(A3) 



The ratio N^__^^+^^_^jN^o^^^^- follows Eq. 
assuming the production fractions are the same for the 
5° and mesons i66]. Table IXiXl lists the estimated 
size of the feed-in background contribution from scmilep- 
tonic decays of b mesons to baryons. 



APPENDIX B: SEMILEPTONIC DECAYS OF 
OTHER b BARYONS 

In addition to the feed-in backgrounds from the 
semileptonic decays of b mesons to baryons, contribu- 
tions are also expected from the semileptonic decays of 
other b baryons. Until recently [1^, [T^I, the A° was the 
only b baryon which had been observed unambiguously. 
Therefore, in order to estimate the number of feed-in 
background events, the production cross section of the 
other b baryons and the branching ratio of the feed-in 
channel must be estimated. The first step in the estima- 
tion is to identify possible contributions to the feed-in. 

Of the lowest lying b baryons, the members of Sf, triplet 
are expected to decay to A^Jtt via the strong interaction 
and contribute to the A° signal. This leaves 5^, 5°, and 



and they are expected to decay predominantly to ^' 
and ilj?. However, by vacuum production of one or more 
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TABLE XIX: Feed-in backgrounds to Aj^ A+ /i from b mesons. All the numbers in parentheses are estimated uncertainties. 
The definition of quantities listed in each column follows Table FVl 

Mode Bi%) (^)^ (SE-)' 

Ag^A+TT- 0.36+[°-^g| TOOO - - 179 ± 19 

K+^i-X - - - - 1237± 97 

Ag^A+^-p^ 7.3 ±(1.4) 0.303 ± 0.004 6.118 1.000 

^A+p^-p^ 0.15 ±(0.15) 0.035 ± 0.002 0.024 0.004 4.3±0.5 

B"^A+nn'i?i, 0.15 ±(0.15) 0.037 ±0.002 0.025 0.004 4.5 ±0.5 



qq pairs, these b baryons can decay into a A^, specifically, 

"6 ^ '^^c ^ P '^M' 

Since tlie decays have a decay topology similar to 
A° — ^ Sc(2455)7r/i~Pp and the fl^ decay topology is sim- 
ilar to that of A° — > A+Tr+Tr"/^"!/^;, a branching fraction 
of 0.39%, 0.39%, and 0.64% (given in Table NTU^ are 
assigned to these three decays, respectively. CDF re- 
ports observing 17 J /ip'Er events in data from an 
integrated luminosity of approximately 1900 pb"-'^ [70[. 
Assuming that the branching fraction for "E.^ — > J /ip'Er 
is similar to that reported by the Particle Data Group 
for A° J/ -0 A 0] , using a generator- level simulation to 
estimate the ratio of acceptances of A'^Kfi~P^ rel- 

ative to that of J/ipE~ (0.2), and scaling by the 

ratio of luminosities for this analysis and the analysis 
(172/1900), each of the Ei, decays is found to contribute 
approximately 0.2% to the signal. For the fl^ decays, a 
similar calculation is performed assuming the same pro- 
duction rate as for the S^. However, because of the larger 
number of particles in the fl^ decay, the acceptance for 
the is an order of magnitude smaller than that of the 
Sfc. For the fl^ , the contribution is 0.03%. These three 
decays are found to contribute < 1% to the signal and 
may be ignored [7l| . 

APPENDIX C: THE bb/cc BACKGROUND 

The bb/cc background refers to the pairing of a A+ 
and a real muon from the decays of two different heavy- 
flavor hadrons produced by the fragmentation of bb or cc 
pairs. In pp collisions, the b and c quarks are primar- 
ily pair-produced via the strong interaction; the single- 
quark production cross-section via the electroweak pro- 
cess, pp W"*" anything — > b anyt hing or c anything, is 
more than 20,000 times smaller |73.[73|. Figure [T71 shows 
the Feynman diagrams up to for the three processes 



which contribute to the bb/cc production [zl, [t^: flavor 
creation, flavor excitation, and gluon splitting. Flavor 
creation, referring to gluon fusion and quark anti-quark 
annihilation, tends to produce bb/cc pairs with an az- 
imuthal angle distribution {A(f>) between the two quarks 
which peaks at 180°. In contrast, the Acf) distribution 
is more evenly distributed for the flavor excitation and 
the low-momentum gluon splitting and peaks at small 
angles for the high-momentum gluon splitting. When 
Acf) is small, daughters of the two heavy-flavor hadrons 
from the fragmentation of bb/cc may appear to come 
from the same decay vertex, as shown in Fig. 1181 If one 
hadron decays semileptonically, and the other hadron de- 
cays into a final state including a A+ — > pK~Tr^ decay, 
the muon from the semileptonic decay together with the 
A+ may be misidentified as the exclusive semileptonic 
signal, AjJ A+/x~j>^. An estimate using pythia has 
shown that this measurement is most sensitive to the 
bb/cc background from high-momentum gluon splitting. 



In the following, the determination of the bb back- 
ground contribution is described. The same procedure 
is followed for the cc background. The ratio N^^i/Nhad is 
given by: 

a,,V (b ^ AtX)V{b ^ ^,-X) e,,_,^+^_^ 
aAo6(Ag^A+^-)e^o^^+^_ 

The CTj,{, is the production cross-section of bb pairs; 
V {b ^ A'^X) and V (b ^ are the probabilities for 

a b and a b quark to fragment into a b hadron and a b 
hadron and to decay to a final state including a A+ and 
a fi~ , respectively. The ^tb^A+fi-x acceptance 
times efficiency for reconstructing the background as the 
A° A'^fi^Df^ signal. The denominator of E g. (ICll) 
can be re-written using the CDF measurement j37| de- 
fined in Eq. (jA3p . the CDF measurement of aB+ [iHl 
assuming ag+ — ago [13, and the world average of 
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S(B0^D+7r-): 
Nti _ 1 

A^had G 

The e^o^yY+^- is determined from a signal simulation 
generated with the bgenerator program [t^, [ill and 
Table IXXl lists the parameters for calculating the denom- 
inator of Eq. (|C2[) . In order to determine the numerator 
of Eq. (jC2p . inclusive 66 events are first generated with 
PYTHIA. The pT of the hard scattering, i.e. the part of 
the interaction with the largest momentum scale, is re- 
quired to be greater than 5 GeV/c. At least one b quark 
must have a px greater than 4 GeV/c and \r]\ less than 
1.5. The value of cr^f,, after applying the kinematic re- 
quirements above, is obtained from pythia, since the sta- 
tus of the a hh m easurements at the Tevatron is still incon- 
clusive [rglTsollsH. Then, the gluon-splitting events are 
filtered and the decays are simulated with evtgen. Only 
events with a /i~ and a AJ which pass the generator-level 
trigger and analysis requirements are considered further. 
Ancestors of the fi^ and the A^t determine whether they 
originate from bb pairs or single b hadrons, and are re- 
trieved by tracing the information from the generator. 
The number of events satisfying these criteria divided 
by the number of generated events gives the product 
r{b^ A+X) V{b^ ^i-X) e(,5^A+M-x- Table IMH lists 
the parameters for the determination of the numerator 
of Eq. (IC2l). 

Table IXXIII lists the estimated ratios, -/Vf,ft/iVhad and 
iVcc/iVhad, based on the values in Tables llodixxll The 
N^^l j,g is found to be only 0.3% of the number of in- 
clusive A+/i~X events. The production of bb and cc 
pairs i n pp collisions has not yet been completely under- 
stood |59l . FtqI [sol . IsH . In order to understand how well 
PYTHIA predicts cr^j and CTcc, an indirect cross-check was 
performed by comparing the differential cross-sections of 
inclusive b hadrons, and in pythia with the CDF 
measurements [H, [tI, [s!] (see Appendix IC 1|) . The dis- 
crepancy between pythia and the data cross-sections 
is generally within 10% for c hadrons and 50% for b 
hadrons, which will be included in the systematic un- 
certainty in Section IVlIII Another cross-check using the 
signed impact parameter distributions of the A+ baryons 



(see Appendix IC 2|) indicates a negligible contribution of 
promptly produced Aji" from cc, which is consistent with 
the above estimate using pythia. 



1. Cross-check of the inclusive b hadron, B , and 
cross-sections 

In order to understand how well pythia predicts 
a^l and Ucc, a cross-check was performed indirectly by 
comparing the differential cross-sections of inclusive b 
hadrons, B+, and D° {da {pp ^ D°X) /dpT etc.) in 
PYTHIA with the CDF measurements [H, [tI, Ull . The 
differential cross-section of D'^ in pythia (see Fig. [T9| , 

for instance, is defined as: da (pp D^X)^'*' /dpx = 
(^cc {N DO /N gen) /'^Pt, where N^o is the number of 
in each pT bin and iVgen is the total number of gener- 
ated cc events. The Apx corresponds to the bin width 
of each px bin, which is the same as that in 
The discrepancy between the pythia and the data cross- 
sections is generally within 10% for c hadrons and 50% 
for b hadrons, which is included in the systematic uncer- 
tainty in Section rVlIII 



2. Cross-check using the signed impact parameter 
distributions 

As an additional cross-check, the signed impact pa- 
rameter distributions (signed do) of the AjJ" baryons with 
respect to the primary vertex, in data and the full simu- 
lation, are also compared. The signed impact parameter 
is defined as do = Q{ro — p), where Q is the charge of the 
particle and ro is the distance between the beam line and 
the center of the helix describing the track in the trans- 
verse plane. The parameter p is the radius of the track 
helix. The full simulation includes the A'^ A'^ii~v^ 
signal and feed-in backgrounds, with relative fractions 
following the estimates in Section fVl Bl An excess of the 
signed do distribution in the region close to zero would 
indicate a significant contribution of the cc background 
in the A+/i~X sample. Figure shows good agreement 
between data and simulation, proving that the promptly 
produced Aji" from cc is a negligible contribution to the 
inclusive semileptonic signals. Figure [20l also shows the 
signed do distributions of £)+ and 
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